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I. 



INTRODUCTION 



A. WHAT IS AROD? 

The U. S. Marine Corps has established the Ground/Air 
Tele Robotic Systems (GATERS) program to develop, test_ and 
field robotic devices for direct use by U. S. Marine 
infantrymen. These robotic systems will be operated under 
the concept of tele robotics which implies that a man, 
stationed remotely from the vehicle, will control the device 
as if he is actually on board. 

The Airborne Remotely Operated Device (AROD) is the 

vehicle being developed through the GATERS program to meet 

the Marine Corps' needs for an airborne, short-range, direct 

support remotely piloted vehicle (RPV) . The GATERS program 

office_[l] states that, from a tactical perspective: 

The AROD is ... designed to allow the front line commander 
to see "over the next hill", out to a distance of two 
kilometers and "around the next corner" in an urban 
environment. . 

This tactical requirement has led to a more technical 
statement of what AROD will be. The GATERS program has 
stated [2 ] : 

The AROD will be a ducted-fan, hovering device capable of 
a forward translation speed of 30 miles per hour. The 
vehicle will carry a fibre optic data link and on board 
cameras to support assigned missions out to a distance of 
5 kilometers. Normal operating radius is 2 kilometers. 

More detailed requirements are listed in Appendix A. 
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B. AROD AND OPTIMAL CONTROL 

The AROD is to be operated by a Marine, untrained and 
unskilled in aerodynamic flight. Therefore, AROD must be 
inherently stable and responsive to simple controls such as 
"go up/down", "rotate right/left", "go forward", "go 
sideways". These instructions will be commanded from a 
joystick, through the control uplink to an on-board control 
system. The control system will interpret the command and 
perform the necessary tasks to maneuver the AROD. The 
design of this controller is the major topic of this work. 
Figure 1.1 is a schematic of the prototype AROD. 




Figure 1.1 Airborne Remotely Operated Device (AROD) 
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Chapter II discusses the AROD system and the equations 
which govern its motion. The equations are complex and the 
system has multiple inputs and multiple outputs. One 
approach to controller design for such a system is to 
simplify the equations describing the system as much as 
possible while retaining as much of the system characteris- 
tics as possible. This is the method developed herein. 

The complex multiple-input, multiple-output AROD system 
is appropriately addressed using optimal control techniques 
with the aid of a digital computer [3]. Optimal control, 
therefore, is the approach taken and developed in this work 
to solve the problem of stability and control of the AROD. 

C. AROD AND RADIO FREQUENCY TRANSMISSIONS 

AROD, as previously described, utilizes a fiber optic 
tether to pass control information and a video signal from 
on board observation devices. Therefore, the need for radio 
frequency (RF) transmissions to or from the vehicle is not 
immediately obvious. In fact, under ideal tactical condi- 
tions, . the RF signature, susceptible to hostile electronic 
countermeasures, would be a most undesirable characteristic 
for a flying vehicle to have. Nevertheless, the fiber optic 
tether will be subject to damage and a back-up system must 
be considered for vehicle control, and for sending back the 
video signal as well. Most importantly, training must be 
conducted prior to any conflict or situation in which AROD 
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is to be used. The fiber optic cable, not immediately 
reusable by the lower echelon ground forces, is expensive 
and, again, susceptible to damage. Therefore, it is 
desirable to provide a low cost method for training opera- 
tors and tactical commanders in the use of the AROD; hence 
the need for an alternative to a fiber optic tether. 

A back-up RF communications link to provide operator 
commands to the AROD has already been designed and imple- 
mented. However, this is an uplink only and provides no 
means for the on board cameras to return their wide 
bandwidth video signal to the ground station. Some 
consideration of this need has been made in the GATERS 
program and a transmission frequency in the 800 to 900 
megahertz range has been tentatively approved. [2] It has 
also been recognized that a video signal via RF might be 
adequate even though it was somewhat degraded relative to 
the high bandwidth video signals possible with fiber optics. 
With these assumptions in mind, Chapter IV deals with the 
design of an on board antenna suitable for transmitting a 
video signal from AROD to the ground station. 

D. SCOPE 

In summary, this work will address two different 
problems associated with AROD: 

(1) Design of a control system for stable flight. 

(2) Design of a radio frequency antenna for the AROD. 



14 



1. Controller Design 



System design of the vehicle has been previously 
undertaken by AROD project engineers. The prototype AROD 
has been built and pertinent information is listed in. 
Appendix A. In Chapter II, the model describing AROD is 
developed. Design of the proposed control system is under- 
taken in Chapter III. 

2 . Antenna Design 

Chapter IV deals exclusively with design and 
analysis of an RF antenna proposed for AROD. The design 
approach is to analyze the requirements and determine a 
basic design which can be refined gradually until an 
acceptable solution is found. The computer program 
Numerical Electromagnetics Code (NEC) , developed by the 
Naval Ocean Systems Center is used extensively to aid in the 
design. No attempt is made to recommend specific hardware 
components in order to implement the antennas. 
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II. AROD DYNAMIC SIMULATION 



The theory behind the equations governing aerodynamic 
flight is well developed. Procedures exist which make the 
modelling of complicated, multiple engine, highly maneuver- 
able aircraft almost a step by step task. Unfortunately, the 
dynamics which characterize AROD differ considerably from 
the standard flight model at a very basic level. Like all 
flying vehicles, AROD has six degrees of freedom of motion: 
three linear displacements and three angular rotations, but 
here the similarities end. This ducted fan device has no 
lift producing surfaces other than the propeller . which 
accounts for 100 percent of its lift and 100 percent of its 
complexity. This chapter traces the development of a 
suitable dynamic model of the AROD, concentrating first on 
the significant gyroscopic contribution of its propeller, 
and then on the complete vehicle. The goal is to provide a 
computer simulation with which an experimental control 
algorithm may be integrated for the purpose of analysis and 
improvement. [4:pp.2-22] 

A. THE GYROSCOPIC AROD 

AROD is a gyroscope. The single propeller rotates about 
the longitudinal vehicle axis to produce a downwash or jet 
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of air through the duct which makes up the AROD body. This 
jet is directed by pitching the body forward, backward, or 
side to side so that the vehicle moves away from the jet in 
the direction desired. To appreciate the control effort 
required to pitch the vehicle in a specific direction 
without causing undesired effects (coupling) in the other 
five degrees of freedom, the equations governing a gyroscope 
need to be understood. 

Figure 2.1 pictures a gyroscope which rotates in the 
same plane as the AROD propeller relative to the coordinate 
system described. This axis orientation will also serve to 
define symbols for later use. For illustration, the body- 
fixed and earth-fixed axes do not coincide in Figure 2.1. 
X,Y,Z and P,Q,R are the body-fixed descriptive geometry for 
linear displacement and angular rotation rate. P,Q, and R 
will also be referred to as roll, pitch, and yaw. X^Y^Z' 
and y?,e,^, are the corresponding earth-fixed displacements 
and rotations. X' can be thought of as altitude. <p,e, and 
are the Euler angles. Velocities will be denoted by U,V, 
and W for rate of displacement in the X,Y, and Z directions, 
and U^V^W' for X^Y^Z'. Pictured also for use in later 
discussion are L,M,N (rolling, pitching, and yawing moments) 
and F x ,F y ,F z (applied force components). 

Angular momentum, h, due to the propeller is defined as 

h r — if w r — ’ x J y "** ^ z (2*1) 
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The rotor lies in the YZ plane of the body coordinates. 
Since the coordinate axes have been chosen so that the rotor 
spins symmetrically about the X-axis, h r is directed only 
along X and h r y = h r z = 0. Equation 2.1 becomes 

h r = I r cj r = ih rx (2.2) 

At this point, a simple calculation will illustrate the 
significant influence the rotor has over the AROD. Refer- 
ence 4 states that in developing the equations of motion for 
aircraft with rotors which exert gyroscopic moments on the 
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body "more often than not, such gyroscopic moments turn out 
to be negligible." In the AROD case, however, I r = .015 

ft 2 -lb(mass) and hover speed of the propeller is w r = 7200 
rpm = 753.98 rad/sec. From Eq. 2.2, h r = 11.3 ft 2 -lb (m)/sec. 
Compared with a nominal total mass of 2.64 lb(m) (85 lb 
weight) , it is clear that the angular momentum imparted by 
the rotor is significant and that gyroscopic effects 'will 
play a large part in modelling the dynamic behavior of AROD. 

Complications arising from these gyro effects can also 
be illustrated. If the AROD body and rotor are considered a 
complete gyroscopic system, then the angular momentum of the 
body must be included in the formulation. Angular momentum 
due to a rotating body is noted as h b = i h bx + jh by + kh bz . 
With the orientation of Figure 2.1, Eq. 2.1 becomes 

h = 1 ( h bx + h r x ) + i h by + k h bz (2.3) 

= i (PI X + I r W r ) + j QI y + k RI Z 

where h is the combined angular momentum of the rotor 
and body. 

Angular moments about the axes of a gyro are discussed in 
the next section and defined by Euler [5:p. 93] as- 

L = h x -h y R+h z Q (2.4) 

M = hy - h z P + h x R 
N = h z -h x Q + hyP 

A simulation of this simple gyro model is given in 
Appendix B. Results are shown in Figure 2.2 for the case of 
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a disturbance about the Y axis (pitch) . Note that the 
actual moments of inertia, listed in Appendix A, were used 






Figure 2.2 Gyroscopic Coupling Effects 

in this first step computer simulation. The time response 
shown here is for the simple case with no rotor accelera- 
tion term (w r = 0) and is the result of a one secon step 
for an angular moment, M, beginning at time = 1 second. The 
amplitude is unimportant and has no effect on the time 
response except to scale the result. The purpose of the 
input was to pitch over the coordinate system by applying a 
torque about the Y axis (Q) . Ideally, a stable, well 
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behaved device would simply pitch over and maintain its new 



attitude. This is not the case for AROD and Figure 2.2 
illustrates the precession and nutation due to a severe 
coupling of the pitch (Q) and yaw(R) axes. These gyroscopic 
effects, then, are the reason a "simple" device such as AROD 
presents such a complex controls problem to the designer. 

B. EQUATIONS OF MOTION 

The basic equations of motion are developed for a 
typical aircraft in References 4 and 5. A combination of 

the two approaches is taken here to arrive at the set of 
equations describing AROD. 

The Euler equations of motion for a free flying rigid 
vehicle with six degrees of freedom are derived from 
Newton's second law. For angular moments applied to the 
coordinate axes, we have 



where h= i h x + j h y + k h z and is the vehicle's total 
angular momentum. 



For linear or translational forces applied along the axes, 



L 



dh 
d t 



(2.5) 





F 




d t 



F, 



y 





( 2 . 6 ) 
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where v = iv + jv + kv and is the total linear transla- 
tional velocity of the center of mass(m) of AROD. 

The dynamic characteristics of AROD are represented by 
Equations 2.5 and 2.6. Other relationships are required to 
solve for each unknown and for each physical force and 
moment acting on the body. The development of a complete 
model is undertaken now for the specific case of AROD using 
actual measurements and experimental data from the proto- 
type. 

1. Angular Momentum 

Equation 2.4 is expanded to become 

h x = i (I x P + I r w r ) + j (I x PR + I r w r R) (2.7) 

- k (I X QP + I r a; r Q) 

hy = — 1 Iy RQ + j Iy Q + k Iy PQ 

h z = i I z QR — j I z PR + k I Z R 
where r . 

w Sh , . 

h = St + u * h 

h = i (I X P + IpWp ) + j Iy Q + k I z R 
w = i P + j Q + k R 

Combining Equations 2 . 4 and 2.7, the AROD angular moment 
equations are : 

L = I X P + (I z - I y ) QR + I r w r (2.8) 

M = I y Q + (I x - I Z )PR + RI r w r 

N = I Z R + (I y - I X )PQ - QIpWp 

The rolling, pitching and yawing moments (L, M, and N) may 
be induced by the aerodynamics of the body and environment 
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or by control inputs through the control vanes and are 
referred to as the "applied moments". They are discussed in 
a later section. Equation 2.8 is next rearranged so that 
the unknown variables (dynamic states of interest) may be 
determined. Roll, pitch, and yaw rates are: 

P - i [(Iy - I* )Q R - Ip“r - L] (2.9) 

- ± x 

Q = i [(I x - I X ) PR - RI r“r - M] 

- L y 

R = T [(I x - Iy) R Q + Qlr^r ~ N] 

X Z 

Already obvious is the mathematical explanation for 
the coupling that was observed in the last section. Every 
combination of pitch, roll, and yaw are present in Equation 
2.9. Pitch, roll, and yaw are each "coupled" to the others. 
Additionally, coupling of the rotor angular acceleration 
with roll will cause the throttle controlling the propeller 
speed to influence these states and further complicate the 
system. 

2 . Force Equations 

The expressions describing the translational forces 
are developed from the body accelerations which are 



II 

X 

> 


i U + j RU - k QU 


(2.10) 


^y = 


-i RV + j V + k PV 




• 

v z = 


i QW - j PW + k W 





where v= iv x + j v y + k v z and is the total body velocity 
vector. 
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Equations 2.6 and 2.10 are combined and rearranged to obtain 
the differential equations containing force terms. 





U = (RV - 


QW) 


+ 


F 

—X 

m 


(2.11) 




1 

§5 

II 

•> 


RU) 


+ 


F 
— y 

m 






W = (QU - 


PV) 


+ 


Fz 

in 





where m is the mass of the vehicle. 

Equations 2.11 contain nonlinear terms as did 
Equation 2.9. These nonlinearities are in the form of 
products of states. The force components F x , F y , F z consist 
of lift from the propeller (F t ) , other aerodynamic influ- 
ences (F a ) , and gravity (g) . F t and F a depend largely on 
experimental data in this development and are taken up 
later. 

Gravity is a constant acceleration causing a 

constant force in the -X’ direction (earth-fixed) . Note 
that small changes in weight due to fuel loss are neglected. 
As the body rotates, the body-fixed axes change relative to 
the earth-fixed axes. Therefore, a coordinate translation 
is needed to relate gravity to each of the X,Y,Z axes 
(g x ,gy,g z ) . The gravity components are given by 

i • g= i g x + j g y + kg z (2.12) 

Any translation from earth-fixed to body-fixed axes can be 
described by a rotation first about the Euler angle <p, then 
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0, and lastly V>* Figure 2.3 illustrates these three 
rotations and the resulting relationships between coordinate 
systems. From Figure 2.3 and Eq. 2 . 12 , the gravity com- 
ponents are 

g x = g(cosV> cos©) (2.13) 

g y = -g(sinV> cos©) 
g z = g(sin0) 

The first-order differential equation for translational 
velocities can now be written as 

U = (RV - QW) + g x + (2.14) 

V = (PW - RU) + g y + 

W = (QU - FV) + g 2 + 

We now have equations describing the six degrees of 
freedom of AROD and the foundation for the model has been 
laid. 

3 . Kinematic Equations of Motion 

The coordinate translation through the Euler angles 
<p r Q,ip resulted in an expression for relating gravity to the 
rotating body. However, the price paid for this necessary 
convenience is the addition of three more variables (<p,0,V») 
to the equations. Therefore, three more equations are 

needed and, if possible, the three variables <p,0,V> should be 
expressed in terms of the existing force or moment vari- 
ables. Fortunately a relationship can be found in the 
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Figure 2.3 Earth to Body-Fixed Coordinate Transformation 
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equation describing the body angular rotations. Angular 
movement of AROD has been expressed in terms of the body- 
fixed angular rates P,Q,R. It can also be described in 
terms of the Euler rates <p , 0, and t/> as 

w = i P + j Q + kR = i 1 <p + j ' © + k • rp ( 2 . 15 ) 



""A coordinate transformation similar to - the' one 
illustrated in Figure 2.3 is possible to relate pitch, roll, 
and yaw with the Euler angle rates. However, some simpli- 
fications can be made with each rotation since <p,0,V> are, in 
fact, the variables of interest. The resulting transforma- 
tion is called the kinematic equations of motion: 



<P 

© 
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sin <p 


cos^tan© 
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-sin<p 


sin<p 


COS<£> 


sin© 


COS© 



(2.16) 



4 . Summary of Equations of Motion 

Nine equations and nine unknowns have been developed 
that completely describe the dynamic behavior of AROD. 
P,Q,R, and U,V,W represent the six degrees of freedom and 
<p,Q,ip, relate the constant force exerted by gravity to the 
body-fixed coordinates. Through the Euler angles, the 
behavior of AROD can also be observed relative to the Earth. 
Specifically, the translational velocities for the earth- 
fixed system U^V^W' can be determined from body-fixed 
velocities U,V,W and <£>/©/0 using the coordinate transforma- 
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tion of Equation 2.13. The applied forces and moments are 
made up of external forces and commanded inputs, and are 
represented by F x , F y , F z , and L, M, N. AROD can be 

pictured as the block diagram in .Figure 2.4. The external 
influences are taken up next. 

C. APPLIED FORCES AND MOMENTS 

The applied forces and moments are a result of the 
commanded inputs and the movement of AROD through the air. 
The commanded inputs can change the rotor speed causing a 
change in the lift generated and a torque on the body. 
Commands can also come in the form of control surface 
displacements which steer the AROD and cause the body to 
rotate about its axes . The movement of AROD through the air 
will create drag. Drag forces can also apply both transla- 
tional forces and angular moments on the body and are a 
result of AROD's aerodynamic characteristics. 

1. Aerodynamic Characteristics 

Two methods exist for determining the aerodynamic 
forces and moments on a vehicle [4:p.lll]: (1) by experiment 

or empirical analysis and (2) by theoretical development. 
The great advantage of (1) is accuracy and the advantage of 
(2) is low cost. The fact that AROD is a ducted fan causes 
a complex interaction of the rotor with the body, and the 
downwash with the control surfaces. Hence, the ducted fan 
compounds the accuracy problem and empirical means for 
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Figure 2.4 Block Diagram of AROD Equations of Motion 
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I 




determining the forces and moments is the most suitable 
approach. Sandia National Laboratories has built a proto- 
type AROD and has collected data based on wind tunnel tests 
of the AROD [6]. This data consists of tabular results 
describing the aerodynamic drag coefficients and downwash 
swirl effect; stability derivatives; and physical measure- 
ments of constants such as weight, moments of inertia, and 
servo gains. This data is listed in Appendix A. The 
tabular data forms a basis from which applied forces and 
moments may be determined once the following flight informa- 
tion is known: 

(1) Weight ratio of actual AROD to 85 lb prototype. The 
data was originally obtained with an 85 lb prototype 
and must be scaled for other weights such as the 
model discussed here which is 76.5 lb. 

(2) Propeller angular speed. Propeller speed is 
determined by setting a hover speed (this model 7200 
rpm) and calculating changes from the throttle 
setting of the engine. 

(3) Propeller angular acceleration. Acceleration is 
based on the difference between speed calculated 
from the throttle and speed at hover. This is an 
experimental approximation arrived at by AROD 
project engineers. 

(4) Translational speed U,V,W. Aerodynamic forces on 

the body are based on the speed and direction of 
AROD (angle of attack) . 

The forces and moments which result are of two 
types: aerodynamic and thrust. The forces and moments due 

to the aerodynamic data are referred to as F ax , F ay , F a z , 
and L a , M a , N a . Those due to thrust and propeller con- 
siderations are F t x , F t y , F t z , and L t , M t , N t . 
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No attempt is made in this work to verify the data. 
The relationships leading to the forces and moments dis- 
cussed in this section are taken verbatim from the informa- 
tion supplied by AROD project engineers, 
a. Angle of Attack 

The tabular data (aerodynamic drag coefficients 
and downwash swirl) result in forces and moments relative to 
the total velocity vector, V t 0 1 , of the body. A relation- 
ship is needed to transform these forces and moments to 
body-fixed axes forces and moments. Angles of attack (a,/?) 
are the standard relationship for this purpose. Angle of 
attack (a) is normally defined as the angle in the XZ plane 
between V t 0 1 and the X axis. Angle of attack (y3) is the 
angle in the XY plane between and the X axis (sometimes 

referred to as the side slip angle). [4] The angle of attack 
used in this model is defined and illustrated in Figure 2.5. 
The resulting relationships are given by 



sina = 



sin 0 = 



W 



V 



tot 

V 



(2.17) 



V 



tot 



where 



v tot= y/ u2 + V 2 + W 2 



b. Aerodynamic Forces 

The forces associated with the tabulated data 
are lift (F t ) , drag (F d ) , and side force (F s ) . They are 
also depicted in Figure 2.5. The transformation from lift, 



31 



X 




a 



x 




p 



Figure 2.5 Angle of Attack Illustration 



drag, and side force is accomplished in like manner to the 
transformation in Figure 2.3 and results in 
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( 2 . 18 ) 



the moments applied to the body axes as a result of AROD's 
movement through air can be derived. These moments are 
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referred to as the aerodynamic angular moments of roll, 



pitch. 


and 


yaw 


(L a , M a , 


N a ) and are 


given as 
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where 

Y y is the yaw moment relative to V t 0 t , 

R r is the roll moment relative to V t 0 t , 

P p is the pitch moment relative to V t 0 t . 

2 . Forces and Moments for Control 

The forces and moments previously discussed are 
based on observation of an AROD prototype and are external 
to the system being designed. A second category is those 
forces and moments which are a direct result of the com- 
manded inputs. As mentioned earlier, the inputs control (1) 
the rotor speed and (2) the displacement of the control 
surfaces . 

a. Moment Due to Ducted Fan Effects 

The AROD in hover (constant rotor speed) closely 
resembles a gyroscope. A gyroscope imparts no torque on its 
axis if it spins with a constant angular rotation. If the 
rotor accelerates (speeds up or slows down) a torque is 
applied to the axis. This torque is accounted for in 
Equation 2.8. However, AROD is also a ducted fan and the 
drag between the rotor tip and the inside body wall creates 
a moment about the X axis (P, roll) . The project engineers 
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for AROD have determined an approximation for this moment 
based on experiment which is given as 

Lr r = .0729 (Ao> r 2 ) (ft-lbs) (2.20) 

where Aw r is the difference between the rotor speed and 
hover, Subscript Tr refers to a dependence on the 
throttle, and .0729 is a constant due to the duct and is 
referred to as in Appendix A. _ 

b. Moments Due to Control Surface Displacement 

Command inputs also include the ability to 
displace the four control vanes into the downwash from the 
duct. As these surfaces are displaced, they impart moments 
about the body axes. The vanes are arranged symmetrically 
as in Figure 2.6. The vanes are displaced by a servo 
mechanism connected directly to the top of each surface. 
Vanes (1) and (3) are operated together as "elevators" and 
impart a moment about the Y-axis (pitch) . Vanes (2) and (4) 
together are the "rudder" and contribute a moment about the 
Z-axis (yaw) . Vanes (1) and (3) displaced in opposite 
directions and (2) and (4) displaced oppositely work as 

"ailerons" to impart a moment about the X-axis (roll) . 

Throughout this work, the term aileron, rudder, and elevator 

will be used as defined here. The actual torque applied by 
each combination of vanes was determined experimentally and 
"constants of effectiveness" were calculated by the AROD 
project engineers. 
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The constants of effectiveness are referred to 
as dimensional stability derivatives when developed theoret- 
ically. They are given the symbols L 8 e , M e e , N r e , for their 
contribution of moments about the roll, pitch, and yaw axes 
due to a displacement by the ailerons, elevators, and 
rudder. They can be thought of as acceleration quantities 
per radian of control surface displacement. The resulting 
relationships, concluding in moments about the three body 
axes, are 
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( 2 . 21 ) 



Lt - Lae^a 
M t — Mg e S e 

N t = N re S r 

where 5 is the displacement of the aileron, elevator 
or rudder into the rotor downwash. L a e , M e e , and N r e 
are given in Appendix A. 

Theoretically, each set of control surfaces will' yield a 
constant of effectiveness about each axis (three constants 
for each aileron, elevator, and rudder) . The symmetry of 
AROD causes negligible cross coupling of control surfaces, 
however, and it is ignored. Note also that the effective- 
ness of each control vane was determined for one rotor speed 
and that the change in downwash associated with the changes 
in rotor speeds is also ignored. 

3 . Servo Equations for Control Surfaces and Throttle 

Commanded inputs sent to AROD include "go for- 
ward/backward", "go side to side", "turn around", and "go 
up/ down " . These commands can be satisfied by pitching the 
AROD forward or sideways (pitch and yaw) for translational 
flight, turning the AROD about its centerline (roll) for 
turning about, and changing the propeller speed for altitude 
control (throttle) . Second order servo motors are used to 
displace each control vane as well as the throttle, 
a. Control Surface Servos 

The servos for each control surface are identi- 
cal. Three equations will describe the operation of the 
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elevators, rudders and ailerons. For each of these equa- 
tions, at least two servos are actually operating at the 
same time on different control vanes. Each servo receives a 
command input (voltage) and results in an angular position 
for the servo. 



where 





= -Hi$ a 


- H 2 5 a 
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H 2 u a 


(2.22) 
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constants , 



u a ,u e , and u r are the servo input. 



b. Throttle Servos 

The servo motor used to open and close the 
throttle is identical to the ones used for the control 
surfaces. An additional first order system describes the 
relationship between the throttle position and the power 
delivered to the rotor (thrust, T r ) . The complete third 
order system is 

s\ = — H, s\ - H 2 5 t + H 2 U t (2.23) 

T r = -T r + .5 (5 t ) 

The power delivered to the rotor, T r , directly relates to 
both the change in rotor speed, Aw r , and the force in the X 
direction. These relationships are 

A w r = T r * Pr ef (2.24) 

AF tx = T r 

where Pr e f is a constant relating the rotor speed to the 
power delivered by the engine. 
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D. ASSEMBLING THE MODEL: PROGRAMMING 

In this section, the equations of motion and the 
force/moment relationships are represented in a block 
diagram of the AROD model developed earlier, and the flow 
chart used to program the model is illustrated. 

1. The Complete Block Diagram 

Figure 2.7 is a block diagram of the complete model. 
For a more detailed explanation of the experimental aero- 
dynamic forces and moments. Appendices A and B should be 
consulted. 

2 . Program Flow Chart 

The logic flow chart developed to program the 
equations and the aerodynamic data is pictured in Figure 
2.8. The IBM Dynamic Simulation Language (DSL/VS) source 
code used is given in Appendix B [7] . 

"V 

E. SIMULATION RESULTS 

Graphical illustration of time response is shown here. 
The four inputs u a ,u e ,u r , and u t are applied alternately as 
step inputs to the AROD model. A caution to the reader 
interpreting these results is appropriate here. The body- 
fixed angle rates P, Q, and R are the dynamic states of 
interest in the control problem. Emphasis should be placed 
on the response of these rates to control input and not on 
the earth-fixed Euler angles <p, ©, and . The Euler angles 
are shown here only as a reference for the reader who wishes 
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Figure 2.7 AROD Model Block Diagram 
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Figure 2.8 AROD Simulation Logic Flow Chart 
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to picture what the vehicle's behavior would be relative to 
the Earth. Confusion can result if one tries to directly 
associate these angles with the roll, pitch, and yaw. By 
way of example (refer to Figure 2.1), the time derivatives 
of <p, 0, and ip are equal to roll(P), pitch(Q), and yaw(R) 
only when <p = 0 and 0 and are small. If the body has 
rotated' 90° ( <p = 90), then Q actually approximates' the 
derivative of ip and R the derivative of 0. 

1. Zero Input Response 

Figure 2.9 is the resulting time response of the 
model when all control inputs are zero. For zero input, all 
control surfaces are parallel with the X-axis and not 
displaced into the downwash. The throttle is set for hover 
speed. Note that the pitch and yaw axes remain stable and 
zero as expected, but the roll axis spins opposite to the 
rotation of the propeller. This precession is a result of 
downwash swirl which was measured and included in the 
tabular aerodynamic data. The propeller downwash is not a 
linear flow of air but has a pattern of vortices which 
results in a force applied to the control surfaces leading 
to a moment about the roll axis. 

2 . Aileron Input 

Figure 2.10(a) is the resulting time response when 
the aileron (5 a ) is displaced and all other control inputs 
are zero as in Figure 2.9. Figure 2.10(a) indicates that 
the ailerons can be positioned to counteract the swirl 
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Figure 2.9 Time Response to Zero Input 
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Figure 2.10 (a) Time Response to Aileron Input 




effect and actually roll the AROD oppositely. These results 
also affirm that the roll axis is uncoupled from the pitch 
and yaw axes as Equation 2.9 suggests. Figure 2.10(b) is an 
illustration intended to show that the roll can be con- 
trolled and set to zero when no roll is desired. In this 
example , an initial displacement of 6 a = -.0499 on the 
aileron control vanes is used and a bias equal' to this 
initial displacement is added to u a . 

3 . Elevator and Rudder Input 

Figure 2.11 is the result when the rudder combina- 
tion of vanes (5 r ) is displaced into the downwash and the 
ailerons are set as in Figure 2.10(b) to minimize roll. 
Note the severe coupling between the pitch and yaw axes, as 
expected. This coupling is a major complication to the 
control process and is dealt with in the next chapter. 
Figure 2.12 is the result when the elevator combination of 
vanes (5 e ) is displaced into the downwash, rudder is reset 
to zero, and the ailerons set as above. The results are 
similar to the rudder observations. In both cases, the roll 
has been affected by the pitch and yaw disturbances, a 
result of the tabular aerodynamic data which changes with a 
and p. The amplitude of the elevator, however, is greater 
than the illustration with the rudder and more clearly shows 
the serious problems that these aerodynamic moments present. 
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Figure 2.10 (b) Time Response: Aileron Displaced 
to Eliminate Swirl Effect 
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Figure 2.11 Time Response to Rudder Displacemnt 
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Figure 2.12 Time Response to Elevator Displacement 



47 



4 . Throttle Input 



Figure 2.13 is the result when the throttle (6 t ) is 
opened to increase the rotor speed and all other inputs are 
zero. The rotor applies a torque about the roll axis and, 
indeed, this influence is seen in the results. The torque 
is directed oppositely to the spin direction (-P) and 
results in a less negative roll. However, the equations' for. 
pitch (Q) and yaw (R) include terms coupled to the propeller 
speed through the aerodynamic moments. The aerodynamic 
moments are determined by angle of attack which is a 
relationship between the body-fixed velocities. As U 
increases, the coupling with pitch and yaw becomes more 
pronounced, the angle of attack changes, and the moments 
determined from the tables change radically. The discon- 
tinuity is best observed by the sharp change in roll, pitch 
and yaw at .75 seconds in Figure 2.13. 

5. Response Without Aerodynamic Moments 

It is insightful to observe the behavior of the 
model with the troublesome aerodynamic moments removed. The 
system behavior without the wind tunnel data supplying the 
moments is much more in keeping with expectations. Figure 
2.14 is the response of this modified system to a rudder 
displacement. Note the similarities of the nutation (Q and 
R) and precession (V>) to the gyroscope illustration of 
Figure 2.2. Figure 2.15 is the response of the same system 
to an opening of engine throttle. This smooth system 
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Figure 2.13 Time Response to Throttle Input 
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Figure 2 . 14 Time Response to Rudder Control 
Aerodynamic Moments Removed 
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Figure 2.15 Time Response to Throttle Input 
Aerodynamic Moments Removed 



51 




response is clearly more in keeping with how we expect the 
actual system to behave. 

The validity of the data was not able to be con- 
firmed nor is it being challenged here. However, these 
results do suggest that the application of the data to the 
system model is a problem that is best addressed by the 
engineers who derive it. 
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III. STABILITY AND CONTROL OF AROD 



The equations of motion describing AROD and the forces 
and moments which act on AROD were developed in detail in 
Chapter II. A computer simulation of these equations 
indicated solid stability in hover, but an inability to 
remain stable once disturbed from the hover state. There- 
fore, a method is needed to maintain stability and drive the 
system to a steady-state condition. Additionally, once a 
steady-state is achieved, it is necessary to steer the AROD 
up/ down, forward/backward, side to side, and make it turn 
about its axis (roll) . Hence, a method is also needed to 
control the AROD and cause it to change from a given initial 
steady-state to another desired steady-state. Stability and 
control are the topics of this chapter. Specifically, an 
algorithm based on optimal control theory is developed in 
the form of the regulator problem to address the stability 
issue. The regulator is then adapted to the control issue 
formulated as the tracking problem of control theory. 
First, however, a summary of the entire stability and 
control issue is offered. Then, a simplification of the 
complete AROD system is made so that standard, well founded 
methods may be used to address the stability and control 
dilemma. 
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A. 



SUMMARY OF THE CONTROL PROBLEM 



The task is to provide an algorithm on which to base a 
control system for AROD. The control system intended for 
AROD relies on inputs provided by an operator to command 
desired outputs. This section discusses the control problem 
based on the inputs, the outputs, and the equations which 
relate them. 

1. Summary of the Complete AROD System 

The result of Chapter II was 18 first order dif- 
ferential equations which are listed in Table 3.1. Also 
listed are the equations for useful quantities which were 
derived in modelling the system. 

Figure 3.1 is a simple, open loop block diagram of 
the AROD system. The block labelled "AROD" represents all 




Figure 3 . 1 AROD Open Loop Model 
the equations listed in Table 3.1. The forces, including 
gravity, and moments acting on the block are external 
influences and are variables in the equations over which the 
designer has no control. The inputs of interest are u a , u e , 
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TABLE 3.1 

AROD DIFFERENTIAL EQUATIONS 

State 

/ 

P = [ (I y - I Z )QR - I r w r + I* + L, ]/I x 

Q = [(I z “ I X )P R - Rlr^r + M3 + M^]/I y 

R = [ (I x - I y )PQ + QI r w r + N a + N t ]/I z 

U = (RV - QW) + F ax /m + F tx /m 

V = (PW - RU) + F ay /m + F ty /m 

W = (QU - PV) + F a2 /m + F tz /m 

cp = P + Qsincp + Rcos<ptan© 

© = Qcos<p - Rsin<p 

i = 0^“ + 

r x sm© cos© 
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S e 
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- H 2 6 e 
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= -HiS* r 
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H 2 u r 


• • 

St 


= -HiS* t 


- H 2 6 t 
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H 2 u t 



u r , and u t . These inputs are voltages applied to servos and 
determine the position of the control vanes in the rotor 
downwash and the throttle position. The outputs of interest 
are <p,Q,ip, and U. If these outputs can be manipulated at 
will, then the rest of the equations, which are derivatives 
of these, can be controlled as well. 
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2 . Optimal Feedback Control 



Control systems may be divided into two basic 
categories: open loop and closed loop systems. Open loop 
control was attempted in Chapter II and is pictured in 
Figure 3.1. Step inputs for pitch (Q) and yaw(R) resulted in 
severe nutation and oscillations. Open loop control is not 
suited to the AROD application. Closed loop (feedback) 
control, then, is the direction taken in this work. 

A feedback control system is one which tends to 
maintain a controlled variable at a fixed steady-state value 
[8]. Therefore, a control system must be able to sense 
variations in the system variables of interest in order to 
apply control to them. For AROD, the outputs are the 
controlled system variables. A scheme to relate the output 
to the input is the feedback loop and is pictured in Figure 
3.2. The feedback matrix, K, must be determined so that a 
given control input always results in the same desired 
output. Note that the reference control vector, r c , does 
not match the commanded input vector, u c , except in steady- 
state. The time when r c # Uj. is called the transient 
response time and is also the time it takes the output to 
reach the desired state. It follows that transient response 
time should be as short as possible. Note also that the 
difference, r c -u c , should be greatest when r c is first 
applied and should decay to zero as the outputs reach their 
desired state. This will result in a smooth transition of 
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Figure 3.2 Output Feedback Model 
the output states from beginning to end. In other words, 
there should be little or no overshoot of the desired state 
by the_ transient states. A good feedback control will 
minimize both the transient response time and the overshoot 
of the desired state. 

A method of feedback control that might minimize 
both the transient response time and the overshoot is 
optimal control. Its name suggests that optimal control may 
provide the optimum solution to the problem of feedback 
control for AROD. However, it is more than semantics which 
drives the choice of controllers. 

First, optimal control lends itself nicely to a 
discrete time solution of the control problem. AROD will 
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employ an on board, digital computer to perform inflight 
stability and control. While a continuous time controller 
can be easily discretized in many cases, design of a sampled 
data controller will simplify the procedure. 

Second, optimal control is ideally suited to state 
variable techniques. A state-space approach is desired 
because of the enormity of the system and the complexity of 
equations. Matrix manipulation utilizing computers will 
prove far less tedious than traditional transfer function 
methods . 

Third, and by far the most important reason to 
select optimal control analysis over classical methods is 
the unique ability of optimal control theory to accommodate 
multiple-input, multiple-output systems such as AROD [3]. 
Four inputs and four outputs have been identified. It will 
be shown in the next section that, although some separation 
of states is possible, severe coupling between other states 
precludes a reduction of the system to single-input, single- 
output systems . 

Fourth and last, optimal control is known to provide 
robust and insensitive solutions to the feedback control 
problem. Assuming that an appropriate performance measure 
is chosen to determine the optimal feedback gain matrix, K, 
the solution can be expected to have a fair degree of 
tolerance to plant model inaccuracies. This section 
develops an approximation of the full AROD model so that 
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troublesome nonlinearities can be dealt with. Therefore, 
some inaccuracies are bound to appear in the new, linear 
model compared with the nonlinear system of Table 3.1. 
Clearly, robustness is not only a desired property of the 
controller, it is an absolute necessity if the controller is 
to be applicable to both the linear and nonlinear model of 
AROD. 

With the advantages of optimal control comes a 
certain drawback: optimal control requires the availability 

of all system states [3:p. 22]. Therefore, the feedback 

scheme pictured in Figure 3.2 must be modified to include 
all system states in the feedback loop. Figure 3.3 satis- 
fies this requirement. However, if the system states 




Figure 3.3 Full State Feedback Model 
are not available from sensor measurements or other means, 
computational observers must also be included to provide an 
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estimate of the hidden states. Observers add an extra 
computational task to the onboard computer, but otherwise 
present no problem to the control designer. 

B. LINEARIZING THE AROD MODEL 

Throughout model development in Chapter II, nonlinear 
relationships were identified. These nonlinearities fall 
into two categories: nonlinear combination of states (e.g. 
Equation 2.8) and discontinuous functions (e.g. table look- 
up of aerodynamic force/moment coefficients) . 

A linear approach to controller design is much preferred 
over a nonlinear analysis. Three major reasons drive this 
preference. First, the theoretical basis for controllers of 
linear systems is well developed and easier to implement. 
Second, in airplane design, the use of linear approximations 
for nonlinear systems has yielded good results for the case 
of steady-state flight and small disturbances therefrom [5]. 
Third, the optimal control solution, K, can often be reduced 
to a single, constant gain schedule for linear, time- 
invariant systems. A constant gain matrix, K, will mean 
that a minimum amount of memory storage will be required of 
the on board computer for control implementation [9]. 

In this section, the equations of Table 3.1 will be 
replaced with linear approximations based on the small- 
disturbance theory, steady-state assumptions, and selected 
physical approximations where necessary. 



60 



1. The Small-Disturbance Theory 



A straight forward discussion of the small distur- 
bance theory is contained in [5:p. 106]. Essentially, if 

AROD is assumed to be in a certain steady-state reference 
condition, then motion of the vehicle consists of small 
deviations from the reference. If this steady-state is 
hover (all translational and angular movement is zero) then 
each of the states in Table 3.1 will be zero or nearly zero. 
This is not. a bad assumption in the case of AROD. The 
flight condition in which AROD will spend most of its time 
is the hover. Only small deviations from the hover will be 
necessary to effect translational or angular movement. For 
example, a translational velocity forward (along the Z' 
axis) of four miles per hour requires only a 15 degree pitch 
of the body. The small-disturbance theory and the steady- 
state hover assumption result in the following simplifica- 
tions : 

(1) Products of states are zero since the product of two 
small values is an extremely small value. 

(2) Derivatives of states are zero since the time 
derivative of a small valued time function is 
extremely small. 

(3) The sine of a state is equal to the state and the 
cosine of a state is equal to one. This is the 
small angle approximation for angles less than 15 
degrees . 

Table 3.2 is the set of equations which result when the 
above assumptions are applied to Table 3.1. Note that much 
of the coupling between states and all of the nonlinear 



61 



TABLE 3.2 



SIMPLIFIED EQUATIONS f SMALL-DISTURBANCE THEORY) 
State 



• 






<p 


= P 




0 


= Q 






= R 




• 

p 


= “(Ir/ix^r + ( L a + L t 


)/Ix 


• 

Q 


= -(I p /I y )w p R + (M a + M 


t)/Iy 


• 

R 


= (I r /I z )u, r Q + (N a + N t 


)/I* 


• 

U 


= g + F ax /m + F tx /m 




V 


= F ay /m + F ty /m 




w 


= F a2 / m + F tz/m 




• ft 


= -H-, 5* - H 2 S a + H 2 u a 




» • 
5e 


= -H-, 6 e - H 2 8 e + H 2 u e 




• • 

5p 


= — H-j 5p H 2 5p H 2 Up 




• • 


= -H 1 5 t — H 2 S t + H 2 u t 





products of states have been eliminated. Note that the 
derivative term of rotor velocity, w r , remains in the 
equation for roll (P) . The term is kept because w r is not a 
state and because it is easily approximated. We wish to 
keep as much similarity between the linear and nonlinear 
models as possible, so physical approximations of 
troublesome quantities are important. 
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2 . Physical Approximations of Forces and Moments 



The force and moment terms in Table 3.2 are made up, 
in part, of functions of aerodynamic coefficients measured 
at discrete AROD flight attitudes in a wind tunnel. The 
data, listed in Appendix A, follow no linear pattern and, 
instead, are discontinuous functions of angle of attack. 
The relationships used to include these data in ' the non- 
linear model were discussed in Chapter II. If a steady- 
state condition is assumed, then some simplification of 
these forces and moments is possible. 

It is not desirable to exclude the aerodynamic data 
completely from the simplified model. An analysis of the 
magnitude of the forces and moments which result reveals 
that they are significant in all flight attitudes. However, 
if a constant hover is assumed then only one set of coeffi- 

cients need be considered. Hover condition is defined by 

- 

AROD project engineers as a=90° and /?=0° [6]. In hover 
condition, the forces and moments become those listed in 
Table 3.3. Note that, for the purpose of Table 3.3, the 
total body-fixed velocity is reduced to U. V and W are 
assumed zero. 

Table 3.3 indicates that each force and moment in 
hover flight can be reduced to a linear function of a state, 
with one exception. Recall from Chapter II, we concluded 
that the rotor downwash created a swirl resulting in a force 
on the control vanes and a moment about the roll 
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TABLE 3.3 

HOVER STATE FORCES AND MOMENTS 
Forces Moments 

Aerodynamic Aerodynamic 



F ax = 0 


L a = U( .036) 

- T r ( .045) - 1.62 


F ay = U ( . 0729 ) 


M a = U ( .81) 


F az = -U (1.08) 


N a = -U (.18) 


Thrust 

F tx = (76.5) + T r (2.125) 


Thrust 

L t = — T r (.5) - 5 a (21.29) 


F ty = 0 


M t = — <5 e (14.51) 


F tz = 0 


N t = -S r (16.68) 



Useful Approximations 
w r = 753.98 rad/sec 

w r = T r (10.472) 



axis. This moment was illustrated in Figure 2.10(a). This 
same constant moment appears in the equation for roll 
moment, L A , in Table 3.3. The compensating factor of .0499, 
determined experimentally in the last chapter, can be 
derived by dividing the constant moment, —1.06, by the 
aileron control vane effectiveness, L a e (-21.29). This 
constant moment will be eliminated in the linear model, but 
will be included in the resulting controller as a constant 
addition to the aileron control, u a . The aerodynamic 
moments about the pitch and yaw axes are zero in hover 
flight. While this poses no problem to the linearizing 
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goal, it means that moments in other flight attitudes have 
not been approximated. Thus the controller resulting from 
the linear model will not consider these external moments. 

3 . The Linear Model Assembled 

The complete, linearized equations are listed in 
Table 3.4. Numerical values for the constant quantities 



TABLE 3.4 

LINEARIZED HOVER STATE EQUATIONS 
State 

i = P 

.0 = Q 

0 = R 

P = U( .0236) - T r (3.738) - S a (21.29) 

Q = U ( .4831) - R ( 6 . 75 ) - 6 e (14.51) 

R = U ( . 108 ) + Q ( 6 . 78 ) - S r (16.68) 

U = T r ( .894) 





(D 

<-o 

* — 
K 

1 

II 


- H 2 S a 


+ 


h 2 




K 


= "Hi 


- H 2 S e 


+ 


h 2 


u e 




= "Hi 


- H 2 6 r 


+ 


h 2 


u r 


• 9 


= -Hi 5* t 


“ h 2 6 t 


+ 


h 2 


u t 



were taken from Appendix A. The next section will assemble 
the equations in a form suitable to determine the optimal 
gain matrix, K. 
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C. FORMULATING THE CONTROL PROBLEM 
1. State Space Representation 

The control system of Figure 3.3 is called a 
regulator. As discussed earlier, the goal of this regulator 
system is to obtain a stable closed loop system which will 
force any arbitrary initial states to zero. The first task, 
then, is to assemble the state equations in a form suitable 
for control analysis. The form desired is given by 

x(t)= Ax (t ) + Bu c (t) (3.1) 

where 

x(t) is the state vector, 

A is the plant matrix, 

B is the control matrix, 

u c (t) is the control vector. 

Note that the A and B matrices do not vary with time. A 
time-invariant system greatly simplifies controller design. 

The 18 states listed in Table 3.4 are the complete 
states of the equations of motion for AROD. However, some 
of the states are unnecessary to the control issue and some 
additional assumptions will aid in determining the best 
control scheme . 

The control problem can be simplified if it is 
recognized that translational velocities over the ground (V' 
and W'. in the earth-fixed horizontal plane) are accomplished 
by pitching the AROD forward/backward or side to side. 
Therefore, the states important to this task are roll(P), 
pitch (Q) , and yaw(R) along with their associated angles. 
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The only translational velocity of interest is the earth- 
fixed vertical velocity or altitude rate. For small pitch 
and yaw angles, then, body-fixed V and W are equivalent to 
earth-fixed V' and W' and will be eliminated from the 

equations. Also, the earth-fixed vertical velocity (al- 
titude rate) will be approximated by the body-fixed veloci- 
ty, U. r " 

A caution about the angles <p,0,ip, is appropriate 
here. The linearized equations conveniently eliminated any 
dependence on the earth-fixed Euler angles. While this may 
be a valid assumption for determining the control solution, 
it changes the meaning of the angle states © and ip. <p 

depends only on the angles © and ip which will be kept small 
in the hover condition, but © and ip both depend on the 

roll (heading) angle, <p, which can range from 0 to 360°. 

Therefore, the linearized angles © b and ip b will be used from 

here on to remind us that these are body-fixed pitch and yaw 
angles and no longer relate the vehicle to the earth-fixed 
coordinate system. 

The complete, linearized system can now be put in 
the form of Eq. 3.1 and is listed in Table 3.5. Figure 3.4 

is the complete signal flow graph for the linear system. 

2 . Stability, Observability, and Controllability 
a. Stability 

A gyroscope is a stable, open loop device. 

Small, bounded disturbances of the roll, pitch, or yaw axes 
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TABLE 3.5 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



STATE-SPACE LINEARIZED SYSTEM 
x(t)= [jP/©b /^b /^e /^t /^a /^e /^r 



A(t) = 

0 0 0 
0 0 0 
0 0‘ 0 
0 0 0 



1 

0 

0 

0 



0 0 
0 1 
0 0 
.894 0 



0 

0 

1 

0 



0 

0 

0 

0 



0 

0 

0 

0 



0 0 .024 0 -3.74 0 
0 0 0 0 -1 0 

0 0 .483 0 0 

0 0 .108 0 
0 0 0 0 



0 -21.29 0 
0 0 0 
0 -6.75 0 -14.51 



0 0 
0 0 
0 0 
0 0 
0 0 



0 

0 

0 

0 

0 



0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 



6.78 0 
0 0 



0 0 0 0 0 
0 0 0 0 0 



0 

0 

0 

0 

0 

0 



0 

0 

0 



0 

0 

0 

0 

0 



0 -157.9 0 
0 0 -157.9 



0 0 0 0 0 0 

0 0 0 0 0 0 _ 

0 0 0 0 0 0 

0-0 0 0 0 0 

0 0 0 0 0 0 

0 .5 0 0 0 0 

0 0 0 0 0 0 

0 -16.68 00000 
0 0 0 1 0 0 0 
0 0 0 1 0 0 

0 0 0 0 1 0 

0 0 0 0 0 1 

0 0 -17.77 000 

000 -17.77 0 0 



0 

0 

0 



0 0 



0 -157.9 0 



0 



0 0 



0 



0 



0 -157.9 0 



0 -17.77 0 

0 0 -17.77 



B(t) = 



000000000000 157.9 000 

000000000000 0 157.9 0 0 

000000000000 0 0 157.9 0 

000000000000 0 0 0 157.9 



Note: t is the transpose operator. 



produce bounded outputs. The outputs may be in the form of 
nutation or oscillations, but it will not grow with time 
once the disturbance is removed. This stability was illus- 
trated for AROD in Figure 2 . 14 where the yaw axis was 
disturbed. The results compared favorably to the response 
of a gyroscope. 

The open loop linearized system may be tested 
for stability also. Since it is linear time-invariant, the 
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Figure 3.4 Signal Flow Graph of the Linear System 




simplified system of Table 3.5 is stable if and only if the 
eigenvalues of the continuous time plant matrix. A, have 
real parts all less than zero [10], Our system fails this 
test for stability since it results in seven out of 16 
eigenvalues equal to zero. Therefore, the linear model has 
lost some of the characteristics of the nonlinear model. 
However, the control algorithm that stabilizes the linear 
system and minimizes the time response and overshoot may 
still be applicable to the AROD. 
b. Observability 

Another relationship to test the applicability 



of the linearized 


plant to 


AROD is 


observability. 


An 


observable 


system 


is one for which 


all states may 


be 


determined 


from observations 


of the 


output. This is 


an 


important 


concept 


since an 


optimal 


controller will 


be 



useless if any system states are not known. To determine 
the observability of our linearized plant, the outputs which 
are available from direct measurement must be known. Refer- 
ence 2 indicates that three single axis rate gyros, one 
vertical rate gyro, a magnetometer (magnetic heading) and a 
barometric altimeter will be available on AROD. Control 
vane position based on a history of the input is also 
assumed to be known. Table 3.6 is a listing of the states 
which will be measured and the device used for each. 

The output of the linear, time-invariant system 

is given by 
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y(t) = Hx(t) 



(3.2) 



where 



y(t) is the output vector, 

H is the measurement matrix 

dimensioned o x n, 
n is the number of states, 

o is the number of measurements. 



TABLE 


3.6 





MEASURED 


STATES 




Measurement Device 


State Obtained 




Magnetometer 


(p (heading) 




Rate Gyro 


U, P, Q, R 




Control History (Servo) 


5 a , 5 e , <$ r , 5 t 






and from Table 3 . 6 



H = 



i 00000000000000 o" 
0001000000000000 
0000100000000000 
0000001000000000 
0000000100000000 
0000000010000000 
0000000001000000 
0000000000100000 
0000000000010000 



Observability of a system due to its measured 
states can be determined from the observability matrix, W 0 , 
given by [11] 



W 0 = [H „ HA „ HA 2 „ HA"' 1 ] (3.3) 

where 

W 0 is a n x mn matrix, 
m is the number of states measured. 



The rank of W 0 (the number of independent columns) must 
equal n if all of the states are observable from output 
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measurements. For the 16 state, linearized AROD system of 
Table 3.5, W 0 has rank 10. Although a control solution 
might still be found assuming all states were available, it 
would not be a realizable controller. Therefore, a dif- 
ferent approach must be taken. 

Figure 3.4 reveals that the roll(P) and the 
throttle (T r ) coupled together and the pitch (Q) and yaw(R) 
states are likewise coupled as expected from the gyroscope 
analogy. However, these two subsystems are coupled only 
through a force/moment relationship dependent on the total 
velocity of the body, V t 0 1 . In the last section, V t 0 1 was 
approximated by U. A further simplification assumes that 
the moment on the pitch and yaw axis due to vertical speed, 
U, is negligible. The complete 16 state linear system can 
now be considered two independent subsystems of eight states 
each as pictured in Figure 3.5. State space representation 
of the two subsystems is listed in Table 3.7. 

Each subsystem has an observability matrix of 
rank eight. Therefore, the states required in the subsystem 
configuration can be obtained from the measurements, 
c. Controllability 

A final test of the linearized model, before any 
effort is expended in actual controller design, is to deter- 
mine whether it is even possible to find a control sequence 
such that a desired final state may be reached. This con- 
cept is known as controllability [10:p. 205]. The control- 
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Figure 3.5 Linear Subsystem Signal Flow Graphs 



73 



TABLE 3.7 



STATE-SPACE LINEAR SUBSYSTEMS 



Roll/Throttle 
XpT (^-) = 

U 
P 
T r 

«. 



Pitch/Yaw 
*qr (t) = 



A> T ” 



0 0 1 0 0 0 

000 .894 0 0 

0 .02 0 -3.7 -21.29 0 



T ~ 



0 

0 

0 

0 

p 

0 

0 



0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 



-1 

0 

0 



0 

0 

0 



.5 

0 

0 



0 

0 

0 

0 

1 

0 



0 -157.9 0 -17.77 



0 

0 

0 

0 

0 

1 

0 



0 -157.9 0 -17.77 



0 

0 



0 

0 



0 

0 



157.91 0 It 
0 157. 9lJ 



©b 


AaR — 


~0 


0 


1 


0 


0 0 


0 


0 ■ 


'Pb 


0 


0 


0 


1 


0 0 


0 


0 


Q 




0 


0 


0 - 


6.75 -14.51 0 


0 


0 


R 




0 


0 


6.78 


0 


0 -16.68 


0 


0 






0 


0 


0 


0 


0 0 


1 


0 






0 


0 


0 


0 


0 0 


0 


1 






0 


0 


0 


0 


-157.9 0 -17.77 


0 


• 




_0 


0 


0 


0 


0 -157.9 


0 - 


17.77 



3a r — 



0 0 0 
0 0 0 



0 

0 



0 

0 



0 

0 



157.9 0 It 

0 157.9 _ 



lability matrix is given by 

W c = [b, AB, A 2 B, ... , A n_1 B] (3.4) 

where 

W c is a nxin matrix, 
i is the number of control inputs. 

The measure of controllability is the rank of W c . For each 
subsystem, the rank is eight. Hence, a control sequence can 
be found for the linearized system. 
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3 . Discretizing the System 



a. Sampling Frequency 

We wish to obtain a discrete time, digital 
control algorithm. Therefore, the system of Table 3.7 must 
be converted to a sampled system. A basis for determining 
the sampling interval, AT, is needed. The sampling theorem, 
based on the Nyquist criterion requires that [11]: 




— 2f M 



(3.5) 



where 

f s is the sampling frequency, 

AT is the sampling period, 

f N is the natural frequency of the system. 

Multiple-input, multiple-output systems do not 
have a single natural frequency. However, if the system is 
considered a set of several single-input, single-output 
systems then a number of different bandwidths are observed. 

The highest bandwidth can be used as the natural frequency 

of the entire system for sampling purposes. In this way, 
four independent systems, with coupled terms set to zero, 
are extracted. The frequency domain transfer functions are 
listed in Table 3.8 along with their bandwidths. A Bode 
analysis yields the frequencies listed. Each of the 
transfer functions also has the describing equation of the 
second order servo in the denominator. Second order 
analysis yields a servo natural frequency of 2 Hz. Since 
the servo frequency is the highest, it will be considered 
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TABLE 3.8 




SIMPLIFIED TRANSFER FUNCTIONS AND 


BANDWIDTHS 








Natural 

Frequency 


Roll 


(s) 


La eH2 


.74 hz 


u a (s) 


s 2 (s 2 +H 1 s+H 2 ) 


Pitch 


0(S) 


Me e H2 


.61 hz 


u e (s) 


s 2 (s 2 +H-] s+H 2 ) 


Yaw 


rl>(s) 


Nr eH2 


. 65 hz 


u r (s) 


s 2 (s 2 +H 1 s+H 2 ) 


U 


U (s) 


-.89H2 


. 03 hz 


u t (s) 


s ( s 2 +H 1 s+H 2 ) ( s+1 ) 



the natural frequency, f N , of the multiple system. The 
sampling frequency is now given as 

f s > 2f M =4 HZ (3.6) 

Equation 3.6 is a lower bound on the sampling 
frequency that will allow the continuous time system to be 
reconstructed once sampled. Selection of the actual 
sampling rate is based on engineering judgment which calls 
for a factor of 10 between the highest frequency component 
and the sampling frequency [12]. A rate of 25 Hz meets this 
criteria and yields a sampling period of AT=.04 seconds, 
b. Analog to Digital Conversion 

The process of conversion is illustrated in 
Figure 3.6. The block labelled D-A is a digital to analog 
signal converter. The block labelled A-D is an analog to 
digital signal converter. The AROD block is the nonlinear 

76 




Figure 3 . 6 Digital System Feedback Scheme 



model which has been linearized here. The next step is to 
combine the D-A, linear AROD model, and the A-D to obtain a 
discrete time representation of the linear, time-invariant 
system. Since the present linear model consists of two 
subsystems, the task is to convert the two systems and 
determine two control systems separately. The discrete 
version of the closed loop system is pictured in Figure 3 . 7 . 
The discrete form of the system Equation 3.1 is given by 

x(k+l)= (J)x(k) + ru c (k) (3.7) 

where 

k = kAT , 
k+l= kAT+AT, 

AT is the sampling period, 

4> is discretized A plant matrix, 

T is discretized B control distribution. 
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U c Ck; 




XCk) 






x(k+l) =0x(k)+ru c (k) 






H . 



=!> yoo 



Figure 3.7 Discretized AROD Feedback 



The discrete plant and control distribution matrices are 
given by [11] 



<j> = e AAT (3.8) 

r = -/^e A * dsB 
o 



where 

e is the natural logarithm operator, 

s is the Laplace operator, 

ds is the derivative with respect to s. 



Equation 3.8 can be solved by using a series 
expansion of the matrix exponential and integrated using a 
numerical scheme on a digital computer. The program called 
ORCONV [13] is used to obtain the discrete time system 
matrices of Table 3.9. 
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TABLE 3.9 



DISCRETIZED SUBSYSTEMS 



^PT ” 



FpT _ 



1 


0 


• 


0400 


-.0030 -.0167 


0 


-.0002 


0 




-.0003 


0 


0 


1 




0 


.0350 0 


.0004 


0 


0 




0 


.00001 


0 


.0009 


1 


-.1466 -.8217 


-.0015 


-.0134 


0 




-.0299 - 


-.00003 


0 


0 




0 


.9608 0 


.0189 


0 


.0003 




0 


.0007 


0 


0 




0 


0 .9010 


0 


.0275 


0 




.0990 


0 


0 


0 




0 


0 0 


.9010 


0 


.0275 




0 


.0990 


0 


0 




0 


0 -4.334 


0 


.4132 


0 




4.334 


0 


0 


0 




0 


0 0 


-4.334 


0 


.4132 




_ 0 


4.334 _ 




<1 


r 


= 










^QR = 




"l 


0 


.0395 - 


.0054 -.0113 


.0012 -. 


0001 . 


00001 




r -.0002 


.00001 


0 


1 


.0054 


.0395 -.0010 - 


.0130 -. 


00001 - 


-.0002 




-.00001 


-.0002 


0 


0 


.9636 - 


.2667 -.5530 


.0879 -. 


0090 . 


0010 




-.0203 


.0016 


0 


0 


.2681 


.9636 -.0768 - 


.6357 -. 


0009 -. 


0104 




-.0014 


-.0234 


0 


0 


0 




0 .9010 


0 


0275 


0 




.0990 


0 


0 


0 


0 




0 0 


.9010 


0 


0275 




0 


.0990 


0 


0 


0 




0 -4.334 


0 


4132 


0 




4.334 


0 


0 


0 


0 




0 0 - 


4.334 


0 


4132 




0 


4.334 



D. THE OPTIMAL REGULATOR 

A review of the linear control problem is appropriate 
here. We now have two discrete-time subsystems, each 
defined by Eq. 3.7. The roll/throttle subsystem is descri- 
bed by the plant matrix (j) p T and control distribution r p T . 
The pitch/yaw subsystem is described by (j) Q R and r Q R These 
matrices were listed in Table 3.9. We also know that all of 
the states in each subsystem are available either by direct 
measurement or computational observers. We wish to find the 
closed loop feedback matrix, K, in Figure 3.7 such that any 



79 



disturbance of any state from the reference hover condition 
(all states zero) will be damped out, returning the states 
to zero in minimum time with minimum overshoot. Since each 
subsystem is controllable, the control law that satisfies 
these requirements is given by [10:p. 771] 

u c *(k) = -Kx(k) (3.9) 

where * denotes optimal in the minimum time, minimum 
overshoot sense. 

Note that K is a constant matrix. We state this goal to 
reduce to a minimum the amount of computer memory required 
to store the feedback gains and, in general, to simplify the 
control algorithm. The validity of this approach is born 
out by experience with linear time-invariant systems [3]. 
The applicability of the resulting controller to the 
nonlinear system, like the previous linearizing assumptions, 
will be determined by the time response of the full AROD 
model once K is found. The following paragraphs will 
provide analytical justification for the steady-state 
approach. 

1. Optimal Control 

Optimal control theory is concerned with minimizing 
a performance criterion chosen by the designer [3:p. 3]. 

Therefore, the primary focus in developing the controller 
should be on a performance measure that quantifies the goals 
of the design effort. These goals are summarized here: 

(1) Minimize the transient response time 
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( 2 ) 

(3) 

(4) 



Minimize the state overshoot 



Determine a constant gain schedule, K. . 

Operate within the physical constraints of the 
system. 

The performance measure typically used in formulating the 
discrete time optimal regulator control solution is 



J= 7X(N)'Cx(N) 



(3.10) 



+ 7 E [x(k)'Q(k)x(k) + u c (k)*R(k)u c (k) ] 

k = 0 



where 



J is the cost function, 

N is the final time step (NAT) , 
k is the time step index, 

Q(k) is the state weighting matrix, 

R(k) is the control cost weighting matrix, 
C is the final state weighting matrix, 
t is the matrix transpose operator. 



The proper choice of C, Q(k) , and R(k) will ensure that the 
control schedule, u c (k) , that minimizes J will meet the 
first two of our design goals and minimize the control 
effort required. 

The solution to Eq. 3.10 for time-invariant systems 
is found from the recurrence relation 

K(N-k) = -[R(N-k) + I^PCk-ljr]- 1 (3.11) 

x r^pfk-ijd) 

and 

P (k) = M (N-k) * P (k-1) M (N-k) + K (N-k) 1 R(N-k) K (N-k) + Q(N-k) 

and 

M (k) = (() - TK (N-k) 



81 



To meet goal (3) , a constant K, it is necessary to 
know that under certain conditions the recurrence relation 
converges to a steady-state for large enough k. The 
conditions sufficient for a steady-state convergence are: 
constant, positive definite weighting matrices Q and R; 
final state weight C=0 ; and u c (k) unconstrained [10]. The 
final time loses its meaning for the steady-state case since 
we assume convergence prior to reaching time-step N. 
Equation 3 . 10 becomes 

* c 

J= £ [xCk^Qxtk) + U c (k) t Ru c (k) ] (3.12) 

k = 0 

where k c is the time step when J converges . 

The recurrence formula can now be written [13] 



and 

and 



K = [R + I^Pr]* 1 r*P<j> (3.13) 

P = M'PM + K'RK + H* QH 
M = <|> - TK 



where the measurement matrix, H, modifies Q so that only 
states of interest are considered. 



The revised performance measure appears to have no 
dependence on the final state. In fact, the solution, K, is 
sometimes referred to as a " subopt imal" gain schedule. 
However, experience shows that for time-invariant systems, 
the suboptimal gain schedule is often quite satisfactory. [9] 
2 . Selecting the Performance Measure 

The weighting matrices Q and R must be selected to 
define the cost function used to determine the optimal 
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gains. For the linear AROD problem, Q and R are both 
diagonal matrices. Each diagonal element of Q corresponds 
to a weight placed on the deviation of a particular state 
from its desired value. Equation 3.13 introduced the idea 
of weighting only certain states utilizing the H matrix. 
The states we are interested in weighting are the angle 
states, the altitude rate, and the control vane displace- 
ments and displacement rates. The angle states and altitude 
rates describe the steady-state of the system. The control 
vanes are constrained in their movement. Therefore, these 
are the states with which we are concerned with controlling 
through the gain schedule. A definition of the measurement 
matrix, H, strictly for the purpose of the cost function, is 



H 



10000000 
01000000 
00001000 
00000 1 0 0 
00000 'o 10 
00000001 



( 3 . 14 ) 



where H applies to both subsystems. 



The non-zero 



elements correspond to <p,U r 6 a ,S t ,8 a ,S t for roll/throttle 

• • 

and © b ,0 b ,5 e ,$ r ,8 e ,8 r for pitch/yaw. 



Now Q is a 6x6 matrix. If all the states of interest are 
weighted equally, Q is given by 



Q = 



iooooo" 
0 1 0 0 0 0 
0 0 1 0 0 0 
0 0 0 1 0 0 
0 0 0 0 1 0 
0 0 0 0 0 1 



( 3 . 15 ) 
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This reduction in cost states will streamline the iterative 
design process of the gain matrix in the next section. 

Each subsystem has two control inputs: u a and u t 

for roll/throttle, and u e , u r for pitch/yaw. A control 

weight matrix which treats each control cost equally is 

«-[;;] ■ 

The cost function resulting in the best gain matrix 
for our purposes may not be the one which weights each state 
and cost equally. Factors such as the units a particular 
state is scaled to and system constraints on other states 
require scaling the diagonal values of the Q and R matrices 
accordingly. It may also be more important to drive a 
particular state more rigidly than others because of the 
dynamics of the system. In the case of AROD, it is most 

■V 

important to drive the angle states and the altitude rate to 
the desired state because the other states, all derivatives 
of these, follow. The time response of the angle states and 
altitude rate is the standard by which the best cost 
function is selected. Therefore, an iterative approach is 
taken, altering Q and R as necessary to meet the following 
specific objectives: settling time of two seconds to less 

than 10% overshoot, and maintaining constrained states 
within limitations. 

The settling time requirement is chosen for two 
reasons. First, the servo response time (second order time 
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constant) is .5 seconds which places a lower limit on the 
response. The equations are coupled which means action by 
one state will force a delayed, but proportional reaction by 
the others. Two seconds for any response on a slow speed 
hovering vehicle such as AROD is considered adequate. 
Second, if the erratic behavior observed in Chapter II is to 
be controlled, then a quick response to disturbances and 
commands is necessary. 

The control vanes are constrained in displacement to 
+.5236 radians (30°) and in velocity to +.8727 radians per 
second. The throttle is constrained to +100 radians and 
+100 radians per second. The iterative design process must 
account for these constraints by weighting the corresponding 
elements of the Q matrix. Appendix A lists constraints. 

Constraints on the servo states imply that the rate 
of change of other states may be restricted as well. A 
constraint on the output states means a limit must be 
applied to the disturbances from which we expect the system 
to recover in the given time. An angular displacement of .1 
radians (6°) is a reasonable disturbance to the system. 
This does not mean that the system is restricted to motions 
within 6° of the hover condition, only that the response 
time for which we design applies up to these limits. 

4 . Iterative Design Process 

We are now ready to form cost functions for each 
subsystem and determine suitable gain matrices. A numerical 
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algorithm for solving Equation 3.13 is discussed in [3:83]. 
A computer program for iterating the solution and converging 
to the steady state is described in [13] and is used to 
determine the gain matrices in this section. 

Once a gain matrix is determined, the time response 
of the linearized, closed loop system can be tested using 
the linear system computer simulation in Appendix C. The 
process repeats until an acceptable time response is 
obtained. 

a. Roll/Throttle Subsystem 

Table 3 . 10 summarizes the results with various 

cost functions. The values for Q are the diagonals of the Q 

• • 

matrix which are the weights on U, 5 a , 5 t , 5 a , and S t , 
respectively. The values for R are the weights on u a and u t 
in that order. 

Iteration (1) is with a Q normalized to account 
for scaling differences between states. Results indicate 
that the normalized cost function is not adequate. A 
shorter response time is desired. 

Iterations (2), (3), (4), and (5) chronicle a 

trial and error design process where each succeeding cost 
function is chosen based on the response of the previous 
ones. A dominant consideration is the coupling between the 
throttle and roll states. Since is undesirable to have the 
throttle used to correct the roll angle, a heavy cost is put 
on maintaining the throttle at the desired state. 
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TABLE 3.10 



ROLL/THRQTTLE SUBSYSTEM COST DESIGN ITERATIONS 



Settling Overshoot 



Cost 

# 


Initial 

Value 


Time 
to U 


% (absolute) 
>o U 


_8a 


Deviation 
5, L 




(1) 


(p =. 1 


4. 


6. 


0% 


(0) 


.02 


0 


.2 


0 


( 2 ) . 


<p =. 1 


1.5 


3. 


12% 


(0) 


.03 ■ 


-.0002 


.35 


.025 


(3) 


p =. 1 


1.5 


2.5 


10% 


(0) 


.03 ■ 


-.004 


.35 


.06 


(4) 


p =. 1 


1 . 


2. 


6% 


(0) 


.07 


-.0007 


1.2 


-.007 


(5) 


p =.l 


1 . 


2. 


3% 


(-.0001) .055 


.0125 


.87 


.3 




U =.l 


2. 


1 . 


(.004) 


5% 


.045 


-3. 


.275 


-50 



Cost Functions 



( 1 ) 



(2) 



Q 

10 

3 

60 

60 

30 

30 



R 

~60" 

60 



Q 

20 

1000 

.1 

.1 

1 

.1 



R 
1 ' 
.1 



(3) 



Q 

50 

2000 

.01 

.01 

3 

.1 



R 
1 ' 
.1 



(4) 



Q 

50 

2000 

.01 

.01 

.2 

.1 



R 
1 ' 
.1 



(5) 



Q 

50 

2000 

.01 

.01 

.4 

.005 



R 



[-os] 



The final iteration results in an acceptable 

time response and satisfies the state constraints. The 
optimal gain matrix, K p T converged in 14 iterations at 
AT=.04 seconds per iteration. Hence, the steady-state 
constant gain schedule is reached after .56 seconds. Since 
the settling time is 1 and 2 seconds, our assumption that a 
constant gain schedule would provide the optimal solution is 
valid. 



b. Pitch/Yaw Subsystem 

Table 3.11 lists results for the iterative 
design of the pitch/yaw subsystem gain matrix. The values 
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for Q are the diagonal elements corresponding to © b , xf) h , S e , 
S r , 5* , and 5 r , respectively. u e and u r are the control 
inputs in R. 







TABLE 3 


.11 














PITCH/YAW SUBSYSTEM COST DESIGN ITERATIONS 












Settling Overshoot 












Cost 


Initial Time % (absolute) 


Deviation 






# 


Value © ib 


© 


±_ 


s n 5 „ 5 


e 


8r. 




(1) 


© =.l 2. 1.5 


10% 


(.013) .035 -.055 


.35 


-.45 


(2) 


© =.l .6 1. 


8% 


(.0125) .08 .045 


.6 


.5 




(3) 


© =.l .75 1.5 


3% 


(.0125) .05 -.065 


.7 ■ 


-.7 








V> =.l 1. 1.5 i 


(-.012) 


2% 


. 065 . 06 


.75 


.6 








Cost Functions 










(1) 


(2) 






(3) 






Q 


R 


Q 




R 




Q 




R 




100 


r x i 


100 




' 1 




150 




1 






100 


L X J 


100 




1 




150 




1 






.01 




.01 








.01 










.01 




.01 








.01 








2 




1 








1 








_ 2 




1 








1 








— 



















The same procedure as roll/throttle for iterat- 
ing through various possible cost functions is used for the 
pitch/yaw subsystem. As in Table 3.10, only selected 
iterations are listed. Insight gained from the roll/ 
throttle iterations was valuable in reducing the steps 
needed for the pitch/yaw design. 

The final gain matrix, K Q R , converged in 13 
iterations, or .52 seconds which is less than the settling 
time of .75 seconds. Therefore, the steady-state gain 
matrix is an optimal solution for this cost function. 
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c. Steady-State Gains 

Table 3.12 lists the optimal steady-state gain 
matrices resulting from the iterative design. A test for 
stability places all eigenvalues of the discrete time closed 
loop subsystems inside the unit circle. Therefore, the gain 
matrix, K, has succeeded in stabilizing the open loop plant. 









TABLE 3.12 














STEADY 


-STATE 


GAIN SCHEDULES 






Roll/Throttle 


















' -1.99 


-.603 


-.717 


.274 


1.75 


.00994 


.140 


. 00016" 




K pt = 


_ -.585 


111. 


-.205 


28.9 


.582 


1.67 


.00962 


.137 




Pitch/Yaw 


















~ -1.61 


-1.61 


-.553 


.320 


1.05 


-.208 


.127 


-.00264“ 




* 

O 

XJ 

II 


_ 1.57 


-1.60 


-.305 


-.566 


.107 


1.23 


.00096 


.130 

























5. Linear System Optimal Regulator Response 

The graphical results discussed here are for the 
optimal gains of Table 3.12 applied to the linear model with 
state constraints. The DSL/VS simulation is listed in 
Appendix C. 

Figures 3.8(a), (b) , (c) , and (d) are the regulator 
time responses for initial conditions of <p( 0) = .l, U(0) = .l, 
© b (0)=.l, and 0 b ( 0 ) = . 1 , respectively. The design goals for 
overshoot and settling time are satisfied for the linear 
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Figure 3.8 (a) Linear System Regulator Time Response 

with Optimal Gains, ¥><0)-.l 
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Figure 3.8 (b) Linear System Regulator Time Response 
with Optimal Gains, U(0)=. 1 
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Figure 3.8 (c) Linear System Regulator Time Response 

with Optimal Gains, 6(0) = . 1 
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Figure 3.8 (d) Linear System Regulator Time Response 
with Optimal Gains, 0(0) = .l 
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system. These results, for which the optimal gains were 
designed, are a reference to which results with the non- 
linear system can be compared. 



E. THE OPTIMAL TRACKER 

The optimal gains derived for the regulator can be 
extended to the tracking problem. The tracking problem 
characterizes the command-response relationship desired to 
steer and maneuver AROD. In theory, several approaches to 
the tracking problem exist. The method used in this work 
utilizes the error state. The control law which applies is 

^ (k) = -K (x(k) - r c (k) } (3.17) 

where r c (k) is the reference (commanded) state vector. 



This control law is pictured in Figure 3.9. The error state 




Figure 3.9 Tracker Control Scheme 
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as it drove the system states to zero in the regulator 
problem. A full state reference vector is required which 
means a commanded value must be specified for every system 
state. The reference vector for the roll/throttle subsystem 
becomes 



a r cp T C^c • U c » c / c « c c c 3 i. 

and for pitch/yaw 

= [e c Ac ,Qc /Rc ,«'r„ ]* (3.19) 

c c c c 



Q R 



where the values of the reference vectors may be time 
varying and are the desired values of the final state. 



Application of Eq. 3.17 to the linear model is described 
in Appendix C. Figure 3.10(a) is the response for desired 
states ^>=©=.1. The response is within the desired limita- 
tions for settling time (2 seconds) and overshoot (10%) . 
Coupling is quite pronounced between the pitch and yaw. 
Figure _3 . 10 (b) is the response to U=.5 and coupling between 
the throttle and roll heading is seen. Finally, Figure 3.10 
(c) is the response to a step input on yaw (V> b = .l). The 
results are within specification. 



F. RESULTS WITH THE NONLINEAR SYSTEM 

The optimal steady-state gains for the linearized AROD 
model are determined in the previous section. The goal in 
designing the control algorithm is to provide fast response 
to disturbances and command inputs, while maintaining 
minimum overshoot of the final state. The requirement for 
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Initial States: <p= 0, U=0 , © b =0, 0 b =0 

Desired States: y?=.l, U=0, © b = .l / 0 b =0 



Figure 3.10 (a) Linear System Tracker Time Response 
with Optimal Gains, Roll and Pitch Command 
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PITCH 
YAW 
ANGLES 
CD 3 RAD) 



ROLL 
ANGLE 
CIO' 3 RAD) 

ROLL_ 

ANGLE 

(5 X K) 3 FI/feEC) 




Initial States: <p=0 , U=0 , © b =0 , t/> b =0 

Desired States: <p=0 , U= . 5 , © b =0 / t/> b =0 



Figure 3.10 (b) Linear System Tracker Time Response 
with Optimal Gains, Altitude Rate Command 
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Initial 


States : 


<P=0 , U=0, 


o 

II 

XI 

0 


V- b =o 


Desired 


States : 


<p= 0, U=0, 


©b=0, 


V> b = « l 



Figure 3.10 (c) Linear System Tracker Time Response 

with Optimal Gains, Yaw Command 
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quickness is driven by the erratic behavior of the full 
model in Chapter II when it is disturbed from the steady- 
state pitch or yaw positions. This section applies the 
linear system controller to the nonlinear model. Results 
are given first for the controller applied to the nonlinear 
model without the pitch and yaw aerodynamic moments and 
second to the full model with all moments and forces. 

1. Nonlinear Model without Aerodynamic Moments 

Figure 3.11(a) is the nonlinear system step response 
for the heading angle, <p. The settling time is 1.5 seconds 
which is the same as for the linear system. However, the 
overshoot is nearly 15% which is larger than the linear 
case. Figure 3.11(b) is the response to a desired altitude 
rate, U', of .1. Note that for the nonlinear model, the 
altitude rate is derived from the Euler angles and it is no 
longer approximated by the body-fixed velocity, U. The 
altitude rate settles in less than 1.5 seconds, but the 
heading angle, disturbed through coupling, settles in 1.75 
seconds with .05 radians overshoot. Still, the system does 
settle quickly and with zero steady-state error. 

Figure 3.12 is the system response for a desired 
state step input on yaw, t/> b , of .1. Severe coupling to <p 
and U' is observed, while pitch, © b , behaves as expected. 
The reality of the linearizing assumptions becomes evident 
here. The results simply indicate that as the body is 
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Figure 3.11 (a) Nonlinear System Tracker Response 

Heading Angle Command, <p 
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PITCH 
YAW ' 
ANGLES 
C D 3 RAO) 



ROLL 
ANGLE 
CIO' 3 RAD) 

ALTITUDE 

RATE 

C5XI0 3 FpfeEC) 




Initial States: 
Desired States: 



<p—0 , U*=0, © b =0 , t/» b =0 
<P=0 , U ' = . 1 , © b =0 ; t/» b =0 



Figure 3.11 (b) Nonlinear System Tracker Response 

Altitude Rate Command, U' 
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PITCH 

YAW 

AN6LES 
CD 3 RAD) 



ROLL 
ANGLE 
CIO' 3 RAD) 

ALTITUDE 

RATE 

C5 X I0 3 FT/5ED 




Initial 


States : 


o 

ll 

b 

o 

ll 

9 - 


©b=0/ 


o 

II 

-O 


Desired 


States : 


o 

II 

D 

o 

II 

9 - 


© b =o. 


^,-.1 



Figure 3.12 Nonlinear System Tracker Response 

Yaw Command 



pitched (or yawed) , some of the altitude holding thrust will 
be lost and redirected to a horizontal velocity. Figure 
3.13 shows that if the pitch or yaw is kept small enough 
(here ij) b = .03 radians), then the altitude rate can also be 
controlled. Figure 3.14 goes a step further and illustrates 
that if a bounded input (pitching the body for 1 second) is 
applied, then a speed across the ground can be realized 
while maintaining the desired heading and altitude rate. 

2 . Nonlinear Model Complete 

Results for the controller applied to the full model 
including the troublesome pitch and yaw aerodynamic moments 
are shown here. Figure 3.15 is the system response to a 
step input on roll (<p=. 1) A steady-state error of -.04 is 
and is due to the downwash swirl moment previously dis- 
cussed. However, coupling to the other states causing the 
instability noted in Chapter II has been eliminated by the 
controller. Figure 3.16 is the result of a step input for 
U' = .5. Coupling of the velocity to the pitch and yaw axes 
causes enough of a disturbance in angle of attack to upset 
the system beyond the constraints of the control vanes and 
the controller is unable to correct it. To illustrate the 
behavior of the data to pitching and yawing, Figure 3.17 
shows the resulting aerodynamic moments, L a , M a , N a , for 
step input for altitude rate. 
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Figure 3.13 Nonlinear System Tracker Response 
Yaw Command for Minimum U’, <p Deviation 
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15CT 



PITCH 

YAW 

ANGLES 
CIO 3 RAD) 



100- 






-SC- 



ROLL 
ANGLE 
CD 3 RAD) 

ALTITUDE 

RATE 

C5 X O 3 FDSEC) 




Initial States: <p=0 , U ' =0 , © b = 0 , ip b =0 
Desired States: ^»=0, U'=0, © b =0, i/> b = .03/0 



Figure 3.14 Nonlinear System Tracker Response 
Yaw Command is a 2 second step 
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1501 
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CIO' 3 RAO) 



ALTITUDE 

RATE 

(5X I0 3 FPSEC) 
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-u 



1 1 1 I 1 

0 12 3 4 5 




Initial States: <p= 0, U‘=0, 0 b =0 , rp b =0 

Desired States: <p=.l, U 1 =0 , © b =0 , V> b =0 



Figure 3.15 Full Nonlinear Model Tracker Response 

Heading Angle Command 
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PITCH 

YAW 

AN6LES 
C I0' 3 RAO) 



ROLL 
ANGLE 
CIO' 3 RAD) 

ALTITUDE 
RATE - 
C5XD 3 F!/$eC) 





Initial States: <p= 0, U 1 =0 , © b =0 , t/> b =0 

Desired States: <p= 0, U'=.l, 0 b =0 , V> b =0 



Figure 3.16 Full Nonlinear Model Tracker Response 

Altitude Rate Command 
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Initial States: <p= 0, U ' =0 , © b =0 , rp b =0 

Desired States: <p= 0, U ' = .1 , © b =0 , 4> b =0 



Figure 3 . 17 Aerodynamic Moments Response to 
Altitude Rate Command 
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3 . Conclusion 



Clearly, the controller is not adequate to deal with 
such large discontinuities as the aerodynamic moments 
present. The control vanes are the limiting factor. 
Therefore, the data or the way it is related to the states 
in the model is not correct. 

'However, we concluded in Chapter II that 'the 
behavior of AROD without these moments closely resembles the 
gyroscopic model that we expected. If the controller is 
judged based on its performance with the nonlinear model 
without the pitch and yaw moments, then the control design 
is valid. The fact that the controller maintains stability 
while directing commanded inputs to steer and maneuver the 
vehicle supports the linear model approach taken in this 
work. 
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IV. RADIO FREQUENCY ANTENNA DESIGN 



A. REQUIREMENTS FOR A RF ANTENNA 

Detailed technical requirements for the AROD video 
downlink have not been specified by the GATERS project 
office as of this writing. Neither have the tactical or 
operational requirements been addressed in detail. However, 
GATERS project quarterly review conferences have undertaken 
discussion of this aspect of AROD and some consensus has 
emerged [2]. It is from these discussions and from the 
stated operational mission of AROD (see Chapter I) that 
assumptions are made on which to base the antenna design. 

Several facts are known about AROD which pertain to RF 
wave propagation. 

(1) Cost and weight are key factors and a simple design 
utilizing inexpensive, readily available materials 
is most appropriate. 

(2) The AROD is intended to be man carried and easily 
readied for launch. Rugged, lightweight 
construction is required. 

(3) A frequency range in the UHF band of 800 to 900 
megahertz (MHz) has been tentatively identified for 
video signal transmission. 

(4) A range of at least five kilometers is required if 
AROD is to be operated at the maximum planned 
distance from the operator. 

(5) An omnidirectional pattern is desired for training 
purposes to ensure that the video signal can be 
received regardless of AROD's orientation. 
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(6) Angular speed of the propeller during hover is about 
7200 rpm (120 Hz) and a three blade, carbon fiber 
composite propeller is to be used. 

(7) The body is made of a carbon fiber composite and the 
frame is composed of aluminum. The forebody is 
fiberglass . 

(8) The AROD contains an on board computer, sensors, and 
a radio receiver for the command uplink. 

This list by no means includes all of the factors which 

contribute to the electromagnetic character of AROD nor does 

it address every possible requirement for a RF antenna on 

AROD. However, these requirements allow a design goal to be 

determined based on the following assumptions: 

(1) The 12 0 Hz disturbance from the rotating propeller 

and the command uplink RF signal will not interfere 
with the 800 to 900 MHz transmission signal. An 

order of magnitude between signals is a common rule 
for non-interfering signals. The propeller (7200 
rpm) and the command uplink (72 MHz) meet this 
constraint. 

(2) For the purpose of electromagnetic analysis, the 
"AROD can be modelled as the cylinder pictured in 

__ Figure 4.1. The cylinder has the outside 

dimensions of the AROD main body surface shown in 
Figure 1.1. 

(3) The carbon fiber body is highly conductive. AROD 

project engineers have stated that the composite 
used has high conductivity characteristics. 

(4) The carbon fiber body will essentially shield out 

electromagnetic interference from the on board 
computer and sensors. 

(5) The antenna will conform to the body of AROD as much 
as possible to prevent breakage while being carried 
and snagging during operation in closed-in areas. 

(6) The antenna will be fixed to facilitate simple 

launch preparation. 



Ill 




3. DESIGN SELECTION 

A "logical point to start the choice of antenna is with 
devices already in use on similar vehicles. AROD is most 
similar to a satellite in that it is free flying yet not 
fast moving and aerodynamic influences are negligible. One 
UHF satellite antenna in common use is the mesh array 
commonly referred to as the umbrella [14]. However, the 
umbrella antenna is neither rugged nor omnidirectional and 
would be a poor choice for AROD. 

Stripline antennas are arrays of slotted strips which 
are surface mounted to a vehicle. The slotted strips are 
cut for specific frequencies and are generally not tuneable. 
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These antennas are lightweight, able to be flush mounted and 
could be made rugged and durable for AROD. Research for 
this report has found only L-band and S-band variations of 
the stripline satellite antenna in use. However, the 
concept of a flush mounted stripline can potentially be 
exploited in the development of the UHF antenna for AROD. 
[15] A simple antenna to employ is a dipole extending 
from the top of AROD or out the bottom. At a transmission 
frequency of 800 MHz, the wavelength is 

A = - f = 14.76 inches (4.1) 

where c = 1.18 x 10 1 0 inches/second. 

A half-wave dipole, for example would extend 7.38 inches 
above or below AROD. This design is not considered here. 
Instead, an antenna which conforms to the AROD body is 
developed to satisfy two of the needs listed above: the 
need for a rugged, fixed antenna and the presence of a 
conducting surface (the body) . 

The need for a rugged, fixed antenna is obvious and will 
not be discussed further. The conductivity of the AROD 
body, however, can possibly enhance the conformal antenna's 
radiation pattern characteristics if the body behaves like a 
reflecting ground plane. This potential and basic antenna's 
which can take advantage of the reflecting body are dis- 
cussed in the next section. 
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C. THEORETICAL DEVELOPMENT 



Two basic designs for the UHF antenna are considered 
here. Referring to Figure 4.1, the first is a vertically 
polarized dipole or dipole array and . the second is a 
horizontally polarized loop antenna. The dipole array and 
the loop antenna are both very well developed in References 
16 and '17. This work does not endeavor to recreate or 
expand on either books ' theoretical development and the 
reader is encouraged to consult these sources if such 
information is desired. Instead, a brief explanation of the 
important characteristics of the dipole and the loop antenna 
is offered as a preface to the practical considerations 
which drive the design process. 

1. Dipole Array 

a. The Half-Wave Dipole 

An omnidirectional far field pattern is the 
goal. For this reason, a dipole which radiates an omni- 
directional pattern is desired. An obvious choice is the 
half-wave dipole. Fed at the center, the current on a half- 
wave dipole will decrease sinusoidally as it nears the end 
points . 

A single half-wave dipole yields the far field 
pattern of Figure 4.2. The horizontal (azimuth) pattern is 
omnidirectional. If the X direction of the dipole cor- 
responds to the X direction of AROD then it is the horizon- 
tal pattern of the dipole which is of greatest interest in 
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Figure 4.2 Far Field Patterns of a Half-Wave Dipole 



obtaining an acceptable solution for AROD. Therefore, the 
emphasis is placed on obtaining an acceptable horizontal far 
field pattern from subsequent antenna designs. 

b. Circular Array of Half-Wave Dipoles 

A dipole which is close to the AROD body will 
have radiation characteristics far different from the free 
space or free standing half-wave dipole. The conducting 
body will, at the very least, block the dipole pattern in a 
direction "behind" or toward the body from the dipole. An 
array of dipoles which conforms to and encircles the body 
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will allow a signal to radiate through 360 degrees in 
azimuth (<p) . Such an array is called a circular or ring 
array [16:p. 352]. A circular array of identical, equally 

spaced half-wave dipoles each excited at its center will 
result in a non-directional far field pattern. Figure 4.3 
is the resulting far field gain pattern of a four dipole 
ring array. Although the pattern is plagued by peaks and 
valleys creating a wide dynamic range in directive gain, it 
is a good starting point for the design. 

The theory and analysis for circular arrays is 
well developed and an array yielding an acceptable omni- 
directional pattern can be derived using the procedure given 
in [16:p. 354]. However, no procedure is found for an array 
with an obstructing body in its center. For this reason, a 
model of the AROD body which approximates its electromag- 
netic characteristics is useful. To such a model, various 
arrays can be added and the resulting far field patterns can 
be analyzed for adequacy. 

2 . Loop Antennas 

The loop antenna is an attractive alternative for 

the AROD application. A loop would naturally conform to the 

body circumference and could be constructed easily. Loop 

antennas are thoroughly discussed in Reference 17. However, 

analytical results are readily derived only for "small" or 

"moderate" sized loops. A moderate size loop is defined as 

2Pi 

b < 2.5 (4.2) 
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where 



Pi = 3.1415962, 
b is the radius of the loop, 

A is the wavelength. 

A circumferential loop on AROD has a radius b > 15 inches. 
At a transmission frequency of 800 MHz, the wavelength A = 
14.76 inches and equation 4.2 becomes: 

b > 6.39 (f = 800 MHz) (4.3) 

Clearly, AROD does not satisfy the requirements for a 
moderate size loop. 

The difficulty with deriving analytical results for 
larger than moderate size loops lies in evaluating the 
Maxwell's Equations. Larger loops such as an antenna around 
the AROD body require a much more extensive mathematical 
development to evaluate the equations than is desired here. 
Therefore, a numerical technique with the aid of a computer 
is best suited for evaluating the radiation pattern of the 
proposed design. 

Figure 4.4 is the far field gain pattern of a single 
loop antenna with a radius b = 15 inches. While this 

pattern is apparently directional and inadequate for the 
AROD application, combined with the AROD body model, a 
better pattern may result. 

3 . Ground Plane Effects 

Perfect ground planes are defined as infinite, 
plane, and perfect conductors. Practical approximations to 
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Figure 4.4 Horizontal Far Field Pattern- Single Loop 
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a perfect ground plane are planes much larger than A 2 and 
good conductors on the order of 10 7 mhos/meter. [16:p. 239] 

Clearly, a cylindrical surface only slightly larger than 2A 
at its greatest dimension can hardly be considered infinite 
and plane. Nevertheless, the fact that the AROD body does 
conduct suggests that some link to ground plane theory can 
be achieved. 

Ground plane theory is the study of reflections. An 
ideal dipole oriented parallel to a perfect ground plane has 
an image in the plane which is directed out from the plane. 
The effect is to double the image radiated outward and 
perpendicular to the plane. In the case of AROD, the image 
is increased outward from the body in the horizontal plane. 
Further, since the body surface is curved, the image will be 
reflected at increasingly obtuse angles as distance from the 
dipole increases. The object of modeling the AROD body is 
to determine if this pattern of reflection can be used to 
enhance the far field patterns of Figure 4.3 and 4.4 and 
obtain an adequate antenna design. 

a. Ground Planes and Phase 

Before proceeding, an important consideration is 
phasing of the signal. The 800 to 900 MHz waveform planned 
for the video signal carries with it a period and phase. 
The phase of this waveform will vary as the distance from 
the antenna increases. Reflection theory stipulates that 
the distance a waveform travels to and from a reflecting 
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surface can be viewed as a straight line distance from the 
original source, in this case the antenna. Thus, a phase 
change is associated with any reflection. This phase change 
can act to complicate or, at 180 degrees, it can cancel the 
original signal. To minimize any such complications, the 
antenna should be positioned flush to the reflecting surface 
or at "a distance resulting in zero phase change in the 
reflected signal. Further work, not considered here, might 
utilize phasing as a means to steer the pattern as desired, 
b. Ground Planes and the Loop Antenna 

One very important result of modelling the body 
as a conductor and approximating it as a ground plane is 
that it will significantly alter the characteristics of the 
loop antenna. In fact, a circumferential loop around the 
cylindrical body might alternatively be considered a line 
source over a ground plane. The "line" is the loop which 
carries a current from its source. Unlike a line source, 
however, this current would be seen on the loop from two 
different directions since the "line" ends are connected 
together. With these important differences in mind, line 
source effects will play a major part in shaping the numer- 
ically generated far field patterns in the next section. 



D. BODY MODELLING AND COMPUTER AIDED DESIGN RESULTS 

The complexity of the body conforming antenna makes use 
of a computer algorithm to model it very practical. 
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Although a simple array of dipoles lends itself to analyti- 
cal solution, a large loop the size of the one proposed does 
not. A computer program written at the Naval Ocean Systems 
Center entitled Numerical Electromagnetics Code (NEC) is 
available on the Naval Postgraduate School ' s mainframe 
computer [18]. This code is used to obtain the far field 
patterns pictured in this work and is very suitable for this 
study . 

The AROD body, although it can be pictured as a simple, 
conducting cylinder (Figure 4.1), is a very complex electro- 
magnetic structure. One approach to simulate such a 
cylinder is to construct a wire grid outline. Another 
method is to utilize a continuous surface model. Both 
techniques are discussed in [16] and both can be modelled 
with NEC. In fact, because of the symmetry of the AROD body 
and the use of symmetry possible in NEC, the continuous 
surface model is somewhat less cumbersome to create. The 
modelling process and the resulting computer code are the 
subject of Appendix D. 

1. Results with Dipole Arrays 

The dipole arrays previously discussed are simulated 
here with the AROD body model. While array effects are 
expected to influence the resulting far field pattern, the 
conducting surface of the AROD will contribute significantly 
as well. All runs utilize a single center feed on a .1 inch 
diameter wire which is divided into 11 equal segments to 
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make a half-wave dipole. The first antenna simulated with 
the body model is the four dipole ring array pictured in 
Figure 4.5. Note that the array radius b = 15.5 inches. 




Figure 4.5 Body Model with 4 Dipole Ring Array 

The half inch beyond the body radius of 15 inches is to 
account for practical considerations of mounting the dipoles 
and to prevent "electrical contact" with the body in the 
computer model. Results of the four dipole simulation are 
shown in Figure 4.6 and should be compared with the free 
standing ring array of Figure 4.3. While the resulting 
pattern is hardly adequate, it does demonstrate the benefi- 
cial effects of the reflecting AROD body. Directivity has 
been reduced. 
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Figure 4.6 Horizontal Far Field - 4 Diple Ring and Body 
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Figure 4.6 indicates that the body has effectively 
duplicated the pattern of each dipole and doubled the 
radiating horizontal pattern. Since the pattern multiplica- 
tion depends largely on the distance between radiating 
sources, an analysis of the spacing is revealing. The 
circumference of the 30 inch diameter AROD is 94.2 inches. 
The antenna radius b = 15.5 inches yields an array radius of 
97.38 inches. The four dipoles, then, are each 24.3 inches 
apart along the circumference. In wavelengths this distance 
is 1.65A. This is not the optimal spacing. The spacing 
sought is pictured in Figure 4.7. If the cylindrical 




Above a Ground Plane 

surface can be ignored and instead thought of as a plane, 
then the interaction of the dipole patterns can be supposed. 
Figure 4.2 showed the horizontal pattern resulting from a 
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single half wave dipole. If the surface reflects, then the 
patterns from two adjacent dipoles will combine. Dipoles a 
half wavelength apart will result in patterns adding exactly 
in phase producing sharp peaks and valleys as seen in Figure 
4.7. Patterns from dipoles a quarter wavelength apart will 
tend to smooth each other. Therefore, spacing close to . 25A 
is sought to reduce the dynamic range. Equal spacing is 
needed to maintain a non-directional array as previously 
discussed. 

A circumference of 97.38 inches and wavelength A = 
14.76 inches yields no whole multiple of .25. Table 4.1 
lists the multiples of quarter wavelength spacing 



TABLE 4.1 

QUARTER WAVELENGTH SPACING 


OF DIPOLES 


Approximate 


Number of 


Actual 


Soacincr 


Dipoles 


Spacina 


A/ 4 


26 


3 . 7 inches 


3A/4 


9 


10.8 inches 


5 A/ 4 


5 


19 . 3 inches 


7A/4 


4 


24.3 inches 


circumference= 97.8 inches 


A/4= 


3.69 inches (800 


Mhz) 



possible on the body. Four dipoles did not yield an 
acceptable pattern. Figure 4.8 is the result of five 
equally spaced half-wave dipoles. The pattern is clearly 
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more smoothed although more complex than the four dipole 
case. The horizontal, vertical and total patterns are shown 
in Figure 4.8, but the total pattern (solid line) is of most 
interest here. In this case the smoothing effects of the 
reflecting body has reduced significantly the sharp drop off 
in directive gain characterized by the four dipole example. 

The case of a nine dipole array deserves some 
practical discussion about feed points. The problem of 
feeding multiple points in phase from one transmitter, as in 
the case of AROD, is easy to overcome on the computer but 
more difficult with actual components. Finite distances 
exist between transmitter and feed point and phasing is 
simplified with fewer feed points. Additionally, power 
losses associated with resistive antennas may become too 
great with multiple loads. Therefore, for the nine dipole 
array, the simulation will consider only three equally 
spaced feed points. 

Figure 4.9 shows the results of this arrangement. 
While not a very omnidirectional pattern, the nine dipole 
case reveals that reducing the number of feed points does 
not increase the directivity and, to the contrary, reduces 
the severity of the peak to valley gain difference of 
dynamic range. With this observation in mind, the five 
dipole antenna which produced a smoother pattern may be 
enhanced with fewer feed points. Figure 4.10(a) is the 
result with only one feed point at <p = 0 on the five dipole 
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Figure 4.8 Horizontal Far Field - 5 Dipoles with Body 




0 




Figure 4.9 Horizontal Far Field - 9 Dipoles with Body 
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array. The far field pattern is greatly improved over the 
five feed case and the dynamic range is greatly reduced. 
Figures 4.10(b) and (c) are vertical cuts over zenith (0) 
for two different azimuth angles. The patterns are non- 
directional and show some promise as a suitable antenna. 

2 . Results with Loop Antennas 

Two variants of the loop antenna with AROD body 
model are studied here. As in the case of the dipole 
arrays, reflections play a major role in improving the far 
field pattern over the free space results of Figure 4.4. 

The loop is a single, circumferential wire pictured 
in Figure 4.11. The first variant has one feed point at <p = 
0. The horizontal pattern which results is shown in Figure 
4.12(a) and differs markedly from the free standing loop 
pattern. While the gain peaks of the free standing pattern 
can be detected in Figure 4.12(a) at <p = 0, 60, 120, 240 and 
300 degrees, the other peaks which give a desirable non- 
directional far field are not a result of the loop theory 
analysis. However, if the antenna is viewed as a line 
source as previously discussed, this pattern is essentially 
an 800 MHz wave propagating along the wire. A peak is seen 
every quarter wavelength with only slight variation at tp = 
180° probably due to opposing currents meeting at this 
point. The currents are not in phase when they meet since 
the distance travelled (half the circumference) is not a 
whole multiple of . 5A. The 97.38 inch circumference fits 
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Figure 4.10 (a) Horizontal Far Field 
5 Dipoles with Body 
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Figure 4.10 (b) Vertical Far Field (<p= 0) 
5 Dipoles with Body 
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Figure 4.10 (c) Vertical Far Field (<p=45° ) 
5 Dipoles with Body 
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Figure 4.11 Single Loop and Body Model 
a 6A line source with a feed at <p = 0. The peaks which are 
seen on the far field pattern of Figure 4.12(a) can be 
viewed as the peaks which result from the line source. 

The far field pattern of the single feed loop is 
further described in Figures 4.12(b) and (c) . For the case 
when (p = 0 , the zenith or vertical pattern resembles that of 
a long (6A) lines source endfire array as described in 
[16:p. 180). Discrepancies between the theoretical endfire 

far field and that pictured here are the result of (1) the 
fact that this antenna is a loop, (2) an actual line source 
lenght of 6.5A, and (3) opposing currents due to the shape 
of the wire as a loop. 
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Figure 4.12 (a) Horizontal Far Field - 1 Loop 1 Feed 
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Figure 4.12 (b) Vertical Far Field (<p= 0) 
One Loop - One Feed 
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Figure 4.12 (c) Vertical Far Field (<?= 45) 
One Loop - One Feed 




The second variant is a single loop with twp feed 
points fed out of phase. The feed points are at <p=0 and 
<p=180 degrees. The horizontal pattern is shown in Figure 
4.13 and has more directivity, than the single feed case. 



E. ANALYSIS & DESIGN CONCLUSIONS 
1. Analysis 

The goal of the design process was to model the AROD 
body and determine some basic antennas that yield omnidirec- 
tional horizontal far field patterns. An iterative approach 
utilizing numerical techniques in the form of NEC was used 
to screen various designs. Only a few of these designs were 
discussed. 

The physical requirement most restrictive tp the 
design was that it conform to the cylindrical AROD body. 
This constraint led to the choice of half-wave dipole arrays 
and a loop. The study did not result in a perfectly 
omnidirectional antenna, but two designs which show poten- 
tial for refinement may have been found. Adequacy in this 
case is judged by the dynamic range that a radio receiver is 
required to have in order to sense an 8 00 MHz signal 
regardless of the orientation, in azimuth, of AROD. In 
other words, the horizontal far field pattern of the AROD 
must have a minimum dynamic excursion between peaks. 
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Figure 4.13 Horizontal Far Field (0=0) 
One Loop - Two Feed 
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2 . 



Conclusions 



The two designs which present the best potential in 
for refinement are the single loop, single feed and the five 
dipole, single feed array. This conclusion is a result of 
their overall, relatively smooth horizontal pattern. Both 
of these antennas are described in terms of horizontal and 
vertical far field patterns in Figures 4.12 and 4.10. 

The most important assumption which was made, 
relative to the resulting patterns, is that the AROD body is 
a good conductor. This assumption is based on the electri- 
cal characteristics of the carbon fiber composite used to 
construct the body. It should be noted that this conducting 
shield is needed to provide the reflecting surface exploited 
by the antenna designs as well as to shield frequency 
interference between the antenna and the electronics 
internal to AROD. 
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V. 



CONCLUSIONS 



A. DYNAMIC MODEL 

A computer model of the ducted fan hovering vehicle, 
AROD, was developed by considering the vehicle as both a 
gyroscope and as a flying vehicle. • Experimental data was 
integrated with the model to simulate the forces and moments 
which act on AROD in flight. Problems with the application 
of some of the data to the computer model were experienced. 
However, the model without the problematic data behaved as 
expected of such a device. 

B. CONTROL 

A linearized approximation of the nonlinear model was 
developed so that optimal control could be used to obtain a 
steady-state gain matrix with the ultimate goal of controll- 
ing the nonlinear computer model. Optimal control proved a 
viable means of determining a constant gain schedule for the 
multiple-input, multiple-output linear system. The result- 
ing controller also found application to the nonlinear model 
with the problematic data removed. 



C . ANTENNA 

Two basic antenna designs were analyzed for use as an 
omnidirectional antenna for a UHF video downlink on the 
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AROD. A five dipole array and a single loop antenna were 
found to have the best potential of the designs which were 
tested using a computer simulation. The body of the AROD 
was found to significantly alter the directivity charac- 
teristics of the antennas. 



D. FURTHER WORK 

In the areas of dynamic modelling and control of AROD: 

(1) The aerodynamic data and its application should be 
verified. 

(2) A training simulator based on the linear model can 
be implemented utilizing the gain schedule derived 
here on small, readily available computers. 

In the area of a UHF antenna: 

(1) Modelling the entire AROD including the forebody and 
testing the designs for variation from these 
results . 

(2) Modelling multiple loop antennas with differing 
phasing in an effort to steer the vertical pattern. 

(3) Analyze the antennas for power requirements as 
weight of batteries and power supplies is critical 
to the AROD. 

(4) Test the designs on the AROD prototype for 
suitability. 
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APPENDIX A 



EXPERIMENTAL DATA 



Empirical data received from Sandia National 
Laboratories is listed here. This data was obtained 
experimentally through measurements and wind tunnel testing. 



1. Weights and Moments of Inertia 



Description 



Name Value (units) 



Weight of Prototype 
Weight of Model 



WtPro 85 (lbs) 

WtMod 76.5 (lbs) 



Moments of Inertia: 



about X-axis 


lx 


1.5246 


(ft 2 lbs) 


about Y-axis 


T y 


1.6767 


(ft 2 lbs) 


about Z-axis 


Iz 


1.6684 


(ft 2 lbs) 


Products of Inertia: (assumed) 


1=1 = 
xy xz 


1=0 

yz 




Moment of Inertia of Propeller 






about X-axis 


Ir 


.015012 (ft^ lbs) 


about Y-axis (assumed) 




0 




about Z-axis (assumed) 


— 


0 




Propeller Weight (derived) 


PropWt 


4 . 50017 


(lbs) 


2 . Hardware Characteristics 








Propeller Efficiency 


Pref 


10.472 


(rad/lb-sec) 


Minimum Thrust of Engine 


Thr min 


35 


(lbs) 


Maximum Thrust of Engine 


Thr max 


115 


(lbs) 


Maximum Throttle Setting 


TA max 


100 


( lbs/sec) 


Maximum Change in Throttle 


RT max 


100 


( lbs/sec* ) 


Propeller Speed at Hover 


HovRat 


7200 


(rpm) 


Propeller Compensation 


RHo 


. 0019 


(ft-lb m ) 


Max Deflection of Vanes 


MaxDef 


.5236 


(rad) 


Max Def Rate of Vanes 


MaxRatDef 


. 87266 


( rad/sec) 
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Description 



Name 



Value (units^ 



Control Surface Effectiveness 
Aileron 
Rudder 
Elevator 



M 

N 



aef f 
eef f 
ref f 



-21.29 (sec 2) 
-14.51 (sec 
-16.68 (sec -2 ) 



3 . Aerodynamic Data 

The data listed in the table below is obtained from wind 
tunnel tests. The manner in which it is applied in the 
simulation (relationship to forces and moments on the AROD 
body) is given in the next section. 

Angle of Attack, AOATOT = >/(a 2 + (3 Z ) 



Name 


50 


55 


60 


65 


70 


75 


80 


90 


CRdel 






mm 


.0807 


.0743 


.0712 


.881 


.085 


Rslope 


. 1514 


. 1533 


. 1673 


. 1977 


. 1834 


. 1764 


.2041 


.20 


Req 


9.293 


8.82 


8.115 


7.848 


7.347 


7.134 


5.011 


0.0 


CPdel 


— 


— 


— 


.3106 


.3081 


.3305 


.3515 


0.0 


Pslope 


.5633 


.5123 


.6974 


.7612 


.7603 


.8189 


.8853 


.90 


Peq 


43.14 


31.00 


29.58 


26.74 


24.55 


20.43 


15.17 


0.0 


CYdel 


— 


— 


— 


. 0024 


-.002 


.0090 


. 0120 


0.0 


Yslope 


. 0609 


. 0663 


. 0215 


. 0060 


-.004 


. 0223 


.0302 


.04 


Yeq 


.555 


.576 


. 006 


-.09 


-1.63 


-1.77 


-1.74 


-1.8 


VanEff 


— 


— 


— 


1.5 


1.5 


1.5 


1.9 


1.0 


Veq 


51.62 


44 . 96 


38.56 


32.94 


27.76 


21.83 


13.31 


0.0 


CLdel 


— 


— 


— 


.2428 


.2391 


.1907 


.1734 . 


1046 


Lslope 


. 7954 


.7183 


.6598 


.5949 


.5901 


.4726 


.4366 


0.0 


Leq 


85. 


85. 


85. 


85. 


85. 


85. 


85. 


85. 


CDdel 


— 


— 


— 


.5754 


.6052 


.5922 


.4966. 


3890 


Dslope 


1.302 


1.313 


1.371 


1.410 


1.493 


1.467 


1.251 


1.2 


Deq 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



CSdel 


— 


— 


— 


.0091 


.0258 


.0331 


.0318 


0.0 


^Sslope 


.0801 


.0801 


.0801 


.0801 


. 0801 


.0801 


.0801 


.0801 


Seq 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



NOTE: ( — ) indicates that no data under this heading was 

received for this angle. 
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4 . Aerodynamic Forces and Moments 

The aerodynamic data listed in paragraph 3, above, 
results in forces moments about the velocity axes, 

The relationships developed by the engineers who obtained 
the data are reproduced here. Some convenient groupings of 
constants and parameters of the system are given here. 



WtRat 



WtMod 

WtPro 



Vdelta = V tQt - V (WtRat) *V eq (AOATOT) 

RHovA = RHo*Pi* (WtRat) * V eq ( AOATOT ) 
w r = angular speed of propeller (rad/sec) 
UseDel = -w r *RHovA*WtRat 
UseEq = (1 + *w r ) *WtRat 

UseSlp = WtRat* (V t t - '(WtRat) *V eq (AOATOT) 
The Forces and Moments which result from the 
aerodynamic data are given here. 

_ EORCES (with respect to 

Due to Body Aerodynamics 

F a l= UseDel *CLdel (AOATOT) + UseSlp*Lslope (AOATOT) 
F ad = UseDel*CDdel (AOATOT) + UseSlp*Dslope (AOATOT) 
F as = UseDel *CSdel (AOATOA) + UseSlp*Sslope 
Due to Thrust Force Aerodynamics 
F-ti= UseEq*Leq 

Total Aerodynamic Forces 

F l= F al + F tl 
F d= F ad 
F s = F as 



147 



MOMENTS (with respect to 

Due to Body Aerodynamics 



R ar = UseDel*CRdel (AOATOT) + UseSlp*Rslope (AOATOT) 
P a p= UseDel*CPdel (AOATOT) + UseSlp*Pslope (AOATOT) 
Y ay = UseDel*CYdel (AOATOA) + UseSlp*Yslope (AOATOT) 
Due to Thrust Force Aerodynamics 
R-tr= UseEq*Req (AOATOT) 

Ptp= UseEq*Peq (AOATOT) 

Y ty = UseEq*Yeq (AOATOT) 

Total Aerodynamic Moments 

R r = R ar + R tr 

R p = p ap + P tp 
Y y = Y ay + Y ty 
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APPENDIX B 



SMJIATTONS SOURCE (DDE 
1. siuple Gyroscope Siimlation 

The following source code is a Dynamic Simulation Language (DSL) 
program to simulate the AROD as a gyroscope. This simulation is 
discussed in Chapter II. 



TITLE BASIC AROO 
* 

CONST I R X = 69 . 5 2 , 

* 

I N l T 

P0= 0 . 

00 = 0 . 

R 0 = 0 . 

MX = 0 . 

MZ = 0 . 



GYRO SIMULATION 

IXX=7063.39, IYY=7768.22, IZZ=7729.58,URX=784. 



DYNAM I C 

IF (TIME .LT. TO) THEN 
MY = 0 
E L S E 

*** STEP INPUT *** 

T 1 = TO + OUR 

MY = AMP*(STEP(TO) • STEP(TI)) 

END IF 

* 

DER l V 



HX 


= I R X * W R X + IXX*P 






HY 


= I Y Y * Q 






HZ 


= I Z Z * R 






PD 


= ( - 1 . *HZ/ I XX )*Q + 


(HY/IXX)*R 


+ MX/ I XX 


QD 


=(HZ/IYY)*P - (HX/IYY)*R + MY/IYY 


RD 


= ( • 1 . * H Y / I ZZ )*P + 


( HX/ I ZZ )*Q 


+ MZ/IZZ 



P = t NTGRL ( P0 , PD ) 

0 = I NTGRL(Q0 , 0 D ) 

R = INTGRL(RO,RD) 

* 

PARAM T 0= 1 , 0 U R = 1 , AMP= 1 0000 

CONTROL FINTIM=5, DELT=.05 
SAVE . 05 , P , 0 , R , MY 
PRINT . 1 , MY , P , 0 , R 
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GRAPH (G0,DE=TEK618) T I M E < N I = 5 , U N = S E C ) M Y ( U N = ' L B - I N ( S Q ) ' ) 

LABEL (GO) INPUT TORQUE 

GRAPH (G1 ,DE=TEK618) T I M E ( N I = 5 , U N = S E C ) P ( U N = ' R A D I A N S / S E C 1 , N I = 3 , . . . 
SC = 2 . , LO= - 2 . ) Q ( L I =3 , N I =3 , SC = 2 . , LO= - 2 . ) R(LI=4,NI=3, ... 

S C = 2 . , L 0 = • 2 . ) 

LABEL G 2 ROLL, PITCH, YAW RESULTING FROM STEP PITCH INPUT' 

END 

STOP 



2. Ocpplete Equations of Motion S.iTTulat-inn 

The following source code is the DSL program used to generate the 
results discussed in Chapter II for the complete AROD model. Notation 
follows that of Appendix ADAIA wherever possible. Data listed in 
Appendix ADAIA is also reproduced here in the form used by DSL. 



TITLE AROD EQUATIONS OF MOTION - FULL MODEL, NO CONTROL - 7 MAY 87 

CONST G R A V= 32.1 74, PI= 3.141 5962 
ARRAY TRANSC9), ABTRAN(9) 



* ★ 
* 

★ 

★ 



* 

★ 

* 

* 

* 

* 

* 

★ 

* 

* 

* 

★ 

★ 

* 

* 

* 

* 

* 

* 

* 

★ 



***************: 



MNEMONICS * 

CR/P/YDEL == ROLL/PITCH/YAW COEFFICIENTS * 

R/P/YSLOPE== ROLL/PITCH/YAW SLOPES OF COEFFICIENT CURVES * 

THRMIN/MAX== THRUST FORCE LIMITS (LBS) * 

R/P/XEQ == EQUIVALENT R 0 L L / P I T C H / Y A W FORCE FACTOR * 

VEQ == EQUIVALENT VELOCITY FACTOR * 

VANEFF == CONTROL SURFACE EFFECTIVENESS * 

CL/S/DDEL == L I FT/SIDEFORCE/DRAG COEFFICIENTS * 

L/S/DSLOPE = = L I FT/S I DE FORCE/DRAG SLOPES OF COEFFICIENT CURVES * 

FRCSLG == CONVERSION FOR I N * * 2 - L B ( F 0 R C E ) TO FT**2-SLUGS * 

WTRAT == RATIO OF AROD WEIGHT TO WIND TUNNEL MODEL WEIGHT * 

p/Q/R == ROLL/PITCH/YAW ANGLE RATES FOR EOM (RAD/SEC) * 

UU/VV/WW == BODY-FIXED VELOCITY COMPONENTS ( X / Y / Z ) ( F T / S E C ) * 

PHI/THT/PSI== EARTH-FIXED EULER ANGLES (RADIANS) * 

ALPHA/BETA/AOAT OT = = ANGLES OF ATTACK (RADIANS) * 

DELE/R/A == ELEVATOR/RUDDER/AILERON/THROTTLE POSITIONS BASED ON * 
IDENTICAL SERVO EQUATIONS (RADIANS) * 

THROT == CHANGE IN THRUST FORCE BASED ON THROTTLE SERVO (LBS)* 

ALT == ALTITUDE BASED ON INTEGRATION OF VERTICAL SPEED (FT)* 

U/V/WERTH == EARTH FIXED VELOCITIES (X/Y/Z) (FT/SEC) * 

GRNSPD/D I ST = = SPEED AND DISTANCE (NO DIRECTION) OVER G R 0 U N D ( F T / S ) * 
THRHOV == THRUST FORCE REQUIRED IN HOVER (LBS) * 

HOVRAT = = ANGULAR RATE OF PROPELLER IN HOVER (RAD/SEC) BASED * 

ON EXPERIMENTAL DATA * 

THR == TOTAL THRUST FORCE (LBS) * 
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== IMPIRICAl DATA FOR VELOCITY CALCULATIONS 
== EQUIVALENT VELOCITY FOR FORCE CALCULATIONS (FT/SEC) * 
== WEIGHTING FACTOR FOR AERO COEFFICIENTS BASED ON * 

VEQ AND AOATOT 

== CHANGE IN ROTOR ANGULAR VELOCITY (RAD/SEC) * 

== ROTOR ANGULAR VELOCITY (DISCRETE) (RAD/SEC) * 



== LIFT/SIDE FORCE/DRAG AERO FORCE COMPONENT (LBS) 

== LIFT/SIDE FORCE/DRAG THRUST FORCE COMPONENT (LBS) 

== AERODYNAMIC FORCES (LBS) * 

= = THRUST INDUCED FORCES (LBS) __ * 

== ROLL/PITCH/YAW AERODYNAMIC MOMENT COMPONENT (FT-LBS)* 
== ROLL/PITCH/YAW THRUST MOMENT COMPONENT (FT-LBS) * 

== ROLL/PITCH/YAW AERODYNAMIC MOMENT (FT-LBS) 

= = ROLL/P I TCH/YAW MOMENT DUE TO THRUST/ROTOR (FT-LBS) 

= = BODY-FIXED GRAVITY COMPONENT ( F T / S E C * * 2 ) 



•»■»************•***•**•*•**•****•******•*■*■***•*■**■»■»■»■»■»■»*■»**■»***■»■»■»***********: 
********** EXPERIMENTALLY OBTAINED DATA *********».„*.», 



★ 


NEUVEQ 


★ 


VDELTA 


* 


R HOVA 


★ 




* 


DE LT I P 


★ 


VT I P 


★ 


USEDEL/S 


* 


DATA 


* 


L/S/DAF 


★ 


L/S/DT F 


★ 


FAX/Y/Z 


* 


FTX/Y/Z 


* 


R/P/YAM 


★ 


R/P/YTM 


★ 


LA/MA/NA 


★ 


LT/MT/NT 


★ 


GX/Y/Z 


★ 





********************************************************************************** 



★ 

★ 


UNITS 








★ 

★ 


* 


MOMENTS OF INERTIA 


I X X , IYY, IZZ, IRX 


LB - FT - SEC 


★ it 2 


•k 


* 


FORCES 


FATX, F A T Y , FATZ, ETC 


LB 




★ 


★ 


MOMENTS 


LAT , MAT , NAT , ETC 


L B - F T 




★ 


★ 


ANGLES 


PS I , THT , PH I 


RAD IANS 




★ 


* 




PSIDEG # THTDEG # PHIDEG 


DEGREES 




* 


★ 


ANGLE RATES 


P,Q,R 


RADIAN / 


SEC 


* 


* 


VE LOC I T I E S 


UU,VV,UW,UERTH, ETC 


FT / SEC 




★ 


★ 




UR 


RADIAN / 


SEC 


* 



* * 
*********************************************************************» 



CONST I XXM = 7063 . 39 , I Y YM = 7768 . 22 , I ZZM= 7729 . 58 , I RXM = 69 . 5 5 2 , WE I G H T = 76.5 

* MEASURED IN**2-LB IN**2-LB IN**2-LB IN**2-LB LB 

* 

CONST P R 0 P W T = 4.500 1 7, PROPEF= 1 0.472, RHO= .00192 

* LB UNITLESS PROPELLER W E I G H T / E F F I C I E N C Y 

CONST T H R M I N = 35., T H R M A X = 115. 

* LB LB MIN/MAX THROTTLE POSITION 

CONST TAMAX= 100., R M A X T = 100. 

* LBS/SEC LB/SEC*SEC MAX CHANGE/RATE OF CHANGE OF THRUST 

CONST MAXD F L= . 5236, R M A X = .87266 

* RADIANS RAD/SEC CONTROL VANE MAX DEFLECTION & RATE 

CONST LAEFF* -21.29, M E E F F = -14.51, NREFF* -16.68 

* ROLL, PITCH, YAW VANE EFFECTIVENESS 

CONST H 1 = 17.77, H2= 157.91, HG1= 219.912, HG2= 24674.126 

* E L E V , RUD, AIL, THROTTLE GAINS/ GYRO GAINS 

* 

CONST WE= 1 . , K = .5 
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* THROTTLE INTEGRATION LOOP GAINS 

CONST A L PM I N = .8727, AL PMAX= 1.5 70796, A 0 A M I N = 50., A 0 A M A X = 90. 

* ANGLE OF ATTACK LIMITS FOR TABLE DATA REFERENCE 



AERODYNAMIC PARAMETERS FROM WIND TUNNEL TESTS 



********** 



* NOTE: 

* 



ALL TABLE DATA IS A FUNCTION OF ANGLE OF ATTACK (TOTAL), 
AOATOT RANGES FROM 50 TO 90 DEGREES. 



AFGEN CRDEL = 



AFGEN RSLOPE = 



****** 


ROLL DATA 


********* 






0807 , 


in 

in 


. 0807, 


60 , 


. 0807 , 


65, 


. 0807 


0743 , 


75 , 


.0712, 


O 

OO 


.0810, 


90, 


. 0850 


15 14, 


55, 


.1533, 


60, 


. 1673 , 


65, 


. 1 977 


1834 , 


75 , 


. 1 764 , 


80 , 


. 2041 , 


90, 


. 2000 



FGEN2 



REQ = 



50, 

70, 

8,2,16, 
9.293,8.5, 
7.347,5 .9, 



50,55,60,65,70,75,80,90, 



8.82,7.1, 

7.134,5.4 



8. 1 15 ,7.2, 
5.011 ,4.67, 



0 , 10 , 
7.848,6 
0 . 0 , 0 . 0 



AFGEN CPDEL = 



AFGEN PSLOPE = 





PITCH DATA 




*******1 




3106, 


55 , .3 1 06 , 


60 


, .3106, 


65 , .3106, 


3081 , 


75 , . 3305 , 


80 


, .3515, 


90 , .0000 


5633 , 


5 5 , . 5 1 23 , 


60 


, . 6974 , 


65 , . 761 2 , 


7603 , 


75 , .8189, 


80 


, . 8853 , 


90 , .9000 



FGEN2 



PE Q = 



8,2,16, 50,55,60,65,70,75,80,90, 0,10, 

34.138,11.69, 31.002,10.45, 29.576,9.72, 26.738,9.44, 

20.427,4.66, 15.168,-0.53, 0.0, -5.92 



24.545,6.78, 



************************* 



YAW DATA 



★ ★★★★★★★★★*★★★** 1 



FGEN2 



YEQ 



******** 

AFGEN CYDEL = 50 , .0024 , 55 , . 0024 , 60 , . 0024 , 65 , . 0024 , 

70 , -. 00 1 5 , 75 , . 0090 ,’ 80 , .01 20 , 90 , .0000 

AFGEN YSLOPE = 50 , . 0609 , 55 , . 0663 , 60 , . 021 5 , 65 , . 0060 , ... 

70 , -.0038, 75 , . 0223 , 80 , . 0302 , 90 , . 0400 

8,2,16, 50,55,60,65,70,75,80,90, 0,10, ... 

0.555 , 0.04 , 0.576, -.5 1 , 0.006, -. 46, - . 0 9 , - . 4 9 , . . . 

- 1 .634 ,- 1.3 , - 1 . 771 , -. 55 , - 1 . 736 ,- 1 . 1 , - 1 . 8, - 1 .2 

CONTROL SURFACE EFFECTIVENESS 

50.1.5, 55,1.5, 60,1.5, 65,1.5, 

70.1.5.75.1.5, 80,1.9, 90,1.0 

VELOCITY COMPENSATING FACTOR 

FGEN2 VEQ = 8,2,16, 50,55,60,65,70,75,80,90, 0,10, ... 

51.616,44.93, 44.955,39.11, 38.564,33.00, 32.935,27.45, 
27.757,2 1 .3 1 , 21 . 829,1 5.06, 1 3.3 14 , 1 1 .65 , 0. 0,0.0 

LIFT DATA 



********** 
AFGEN VANEFF 



AFGEN 


CLDEL 


= 


50, 


.2428, 


55 , . 2428, 60 , . 2428, 


65 , . 2428 








70, 


.2391 , 


75 ,. 1 907, 80 , . 1 734 , 


90 , . 1 046 


AFGEN 


LSLOPE 


= 


50, 


.7954, 


55 , . 71 83 , 60 ,. 6598, 


65 , .5949 








70, 


.5901, 


75, .4726, 80, .4366, 


90, . 0000 


CONST 


LEG 


- 


85 . 


******* 


DRAG DATA ********** 




AFGEN 


CDDEL 


= 


50, 


. 5754 , 


55 , . 5 754 , 60 , .5754 , 


65 , . 5754 








70, 


.6052, 


75 , . 5922, 80, . 4966, 


90 , .3890 


AFGEN 


DSLOPE 


= 


50, 


1.3017, 


55,1.3131, 60,1.3711 


, 65,1. 








70, 


1 .4933 , 


75 , 1 .4674 , 80, 1 . 2507 


, 90,1. 


CONST 


DEQ 


= 


0 . 0 
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***************1*********** SIDE FORCE DATA * * 
AFGEN CSDEL = 50, . 009 1, 5 5, .0091, 60, .0091, 

70 , .0258, 75 , .0331, 80, . 031 8, 

CONST SSLOPE = 0.0801 

CONST SEQ = 0.0 

*********************************************** 

INITIAL 

* 

*** CONVERSIONS 

***** FORCE TO SLUGS AND I N * * 2 TO F T * * 2 
FRCSLG= 1 . / ( GRAV* 1 44 . ) 

IXX= IXXM*FRCSLG 
I Y Y = IYYM*FRCSLG 
IZZ= IZZM*FRCSLG 
I R X = IRXM*FRCSLG 



★ ★ ★ 


MOMENTS 


OF 


INERT I A 




I YM I Z = 


I Y Y 


- I ZZ 




I ZM I X = 


I ZZ 


- I XX 




I XM I Y = 


I XX 


- I Y Y 



M A S S = WE I GHT/GRAV 
WT R A T = WEIGHT/85. 

********************** INITIAL CONDITIONS ************************** 
PO = 0 . 

QO = 0. 

RO = 0 . 

UO = 0 . 

VO = 0 . 

WO = 0 . 

PH I 0= 0 . 

T H T 0 = 0 . 

P SJ^O = 0 
D E L E D 0 = 

DELEO = 

D E L RD 0 = 

DELRO = 

D E LAD 0 = 

* DELAO = 

D L T R D 0 = 

DLTRO = 

T H R 0 T 0 = 

ALTO = 

G R N S P 0 = 

T H R H 0 V = 

H 0 V R A T = 

************ 

MACRO ANGOUT 
AN S I GN = 

A N GOU T = 

A N GOU T = 

ENDMAC 

*****•••• DEFINE MACRO FOR USE IN SERVO EQUATION LIMITING 



0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

WEIGHT 

(THRHOV - PROPWT ) *PROPE F 
DEFINE MACRO TO LIMIT ANGLE 
= ANL I M < ANG I N , AN GM I N , AN GMAX ) 

S I GN < 1 . , ANG I N ) 

LIMIT(ANGMIN,ANGMAX,ANSIGN*ANGIN) 

ANSIGN*ANGOUT 



*************************** 



*********************** 

65 , .0091 , ... 

90 , . 0000 
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MACRO OU T = SERVOCOUTO , OUTD , L I M ) 

01) T = INTGRL(OUTO,OUTD) 

OUT = LIMIT(-LIM,LIM , OUT ) 

ENDMAC 

♦A******************************************************************** 

DYNAMIC 



COORDINATE TRANSLATION COMPONENTS 
***** BODY-FIXED TO EARTH FIXED (EULER ANGLES) 

* ROW 1 

TRANS( 1 ) = COS < PS I ) *COS < T HT ) 

TRANS< 2) = COS(PSI)*SIN(THT)*SIN(PHI) 

TRANS <3 ) = COS(PSI )*SIN(THT)*COS(PHI ) 

* ROW 2 

TRANS(4)= S I N < PS I ) *COS < T HT ) 

TRANS(5)= SIN(PSI )*SIN(THT)*SIN(PHI ) 
TRANS(6)= SIN(PSI )*SIN(THT)*COS(PHI ) 

* ROW 3 

TRANS( 7) = - S I N < T HT ) 

TRANS(8)= COS < T H T ) *S I N < P H I ) 

TRANS(9)= COS < T H T ) *COS < P H I ) 



S I N ( P S I ) * CO S ( P H I ) 
S I N ( P S I )*SIN(PHI) 



COS(PSI )*COS(PHI ) 
S I N ( P H I )*COS(PSI ) 



******* VELOCITIES 
VTOT = SORT (UU*UU + VV*VV + WW*WW) 

U E R T H = UU*TRANS( 1 ) + VV*TRANS(2) + WW*TRANS(3) 
VERTH= UU*TRANS<4) + VV*TRANS(5) + WW*TRANS(6) 
WE R T H = UU*TRANS(7) + VV*TRANS<8) + WW*TRANS(9) 
GR N SPD = SORT (VERTH*VERTH + WERTH*WERTH) 



* 

* 

* 



******* ANGLE of attack 



IF(VTOT.LE..5) THEN 
ALPHA= 1 .570796 
BETA 



0 . 



ELSE 



ALPHA= AS I N < WW/VTOT ) 

BETA = AS I N < VV/VTOT ) 

END I F 

ALPHA = ANLIM(ALPHA,ALPMIN .ALPMAX) 

BETA = AN L I M < BETA , 0 1 . 570796 ) 

AOATOT= RADEG*SQRT (ALPHA*ALPHA + BETA*BETA) 
AOATOTs LIMIT(AOAMIN,AOAMAX,AOATOT) 



* * * 
* 



ANGLE OF ATTACK COORDINATE TRANSFORMATION 
BODY-FIXED TO VELOCITY 
ROW 1 

ABTRAN(1)= -COS(ALPHA)*COS(BETA) 

ABTRAN(2)= SIN(BETA) 

ABTRAN(3)= S I N < A L P H A ) * C 0 S ( B E T A ) 

ROW 2 

ABTRAN(4)= CO S < A L P H A ) * S I N < B E T A ) 

ABTRAN(5)= COS(BETA) 

ABTRAN(6)= - S I N < A L P H A ) * S I N ( B E T A ) 
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* ROW 3 

ABTRAN(7)= - S I N ( ALP H A ) 
ABTRAN<8)= 0. 
A8TRAN(9)= -COS(ALPHA) 



******************** FORCE COMPONENTS ****** 
THR = L I M I T ( T HRM I N , THRMAX, T H R H 0 V + T H R 0 T ) 
N E W VE Q = FGEN2< VEQ , AOATOT , 0 . ) 

VD E L T A = VTOT- S Q R T ( W T R A T ) * N E W V E Q 
RHOVA = RHO*P I*SQRT (WTRAT )*NEUVEQ 

* 

BELT IP* (THR - WEIGHT)*PROPEF 
VTIP = HOVRAT + DELTIP 



U S E D E L = -BELT I P*RHOVA*WTRAT 
USESLP= VDELTA*WTRAT 

USEEQ = (1. + ( 2 . / VT i P )*DE LT I P ) * WTRAT 

*** FORCES 

LAF = USEDEL*AFGEN(CLBEL , AOATOT) + U S E S L P * A F G E N ( L S L 0 P E , A 0 A T 0 T ) 

SAF = USEDE L* A F GEN ( CSD E L , AOATOT ) ♦ U S E S L P * S S L 0 P E 

OAF = USED E L *A FGEN ( CODE L , AOATOT ) + U S E S L P * A F G E N ( D S L 0 P E , A 0 A T 0 T ) 

LTF = USEEQ*LEQ 



FAX 


= 


A B T R A N ( 1 ) * D A F 


+ ABTRAN ( 2 ) * S A F 




FAY 


= 


ABTRAN(4)*DAF 


+ ABTRAN(5)*SAF 




FAZ 


= 


ABTRAN(7)*DAF 


+ A8TRAN(8)*SAF 




FTX 


= 


ABTRAN ( 3 ) * L T F 






FT Y 


= 


ABTRAN(6)*LTF 






FT Z 


= 


ABTRAN(9)*LTF 






FAT X 


= 


FAX + FTX 






FAT Y 


s 


FAY + FTY 






F-AJLZ 


= 


FAZ + FTZ 






MOMENTS 











ABTRAN(3)*LAF 
ABTRAN(6)*LAF 
A B T R A N ( 9 ) * L A F 



RAM = USEDE L*AF GE N ( CRD E L , AOAT OT ) + U S E S L P * A F G E N ( R S L 0 P E , A 0 A T 0 T ) 



* 



+ USEEQ* FGEN2( REQ , AOATOT , 0 . ) 

PAM = USEDE L* A F GE N ( CPDE L , AOA TOT ) ♦ U S E S L P * A F G E N ( P S L 0 P E , AOA T OT ) . 

+ USE EQ* FGEN2 ( PEQ , AOATOT , 0 . ) 

YAM = USEDEL*A FGEN ( CYDE L .AOATOT ) + U S E S L P * A F G E N ( Y S L 0 P E , A 0 A T 0 T ) . 

+ USEEQ* FGEN2 ( YEQ , AOATOT , 0 . ) 

LA = ABTRAN ( 1 ) *RAM + A B T R A N ( 2 ) * P A M + A B T R A N ( 3 ) * Y A M 
MA = ABTRAN ( A ) *RAM + A B T R A N ( 5 ) * P A M +ABTRAN (6)*YAM 
NA = ABTRAN( 7)*RAM + A B T R A N ( 8 ) * P A M + A B T R A N ( 9 ) * Y A M 
LT = • OE LT I P*DE LT I P* FRCSLG* . 0729 ♦ I X X * L A E F F * D E L A 
MT = MEEFF*IYY*DELE 
NT = NREF F* I ZZ*DELR 
LAT » LA + L T 
MAT » MA + MT 
NAT = NA ♦ NT 
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★★*★★*★★****★★***★★★★★*★★ GRAVITY COMPONENTS 
GX = - GRAV*TRANS( 1 ) 

GY = -GRAV*( -TRANS(4) ) 

GZ = -GRAV*( -TRANS(7)) 

* i 



*************** 



SAMPLE 

**** INSERT CONTROLLER HERE 
★ 

★ 

DER I VAT I VE 

★ ^ ■ 



MOMENT DER I VAT IVES 



PD 


= 


( I YM I Z 


* R 


* Q - 


I RX 


*DELT I P + 


LAT )/ I XX 


QD 


= 


( I ZM I X 


*P 


*R - 


I RX 


* V T I P * R + 


MAT ) / I Y Y 


RD 


= 


( I XM I Y 


*P 


* Q + 


I RX 


* V T I P * Q + 


NAT )/ I ZZ 




r*********** 


FORCE 


DER I VAT IVES 




UD 


= 


( VV*R 


- 


W W * Q ) 


+ 


F ATX/MASS 


+ 


GX 


VD 


= 


(UU*P 


- 


U U * R ) 


+ 


FAT Y/MASS 


+ 


GY 


UD 


= 


( UU*Q 


- 


VV* P ) 


+ 


FATZ/MASS 


+ 


GZ 






ANGLE DERIVATIVES ************** 

P H I D = P + (Q*SIN(PHI) + R*COS(PH I ) )*TAN(THT ) 

T H T D = Q*COS(PHI ) - R*SIN(PHI) 

P S I D = (Q*SIN(PHI ) + R*COS ( P H I ) ) /COS ( T HT ) 

SERVO DERIVATIVES ******************** 



ELEVATOR , 


RUDDER , 


AILERON, THROTTLE 


D E L E D D = 


- H 1 *DELED 


- H 2 * D E L E 


+ 


H 2 


* U E 


D E L R D D = 


- H 1 *DELRD 


- H 2 * D E L R 


+ 


H 2 


* U R 


DE L ADD = 


- H 1 *D E LAD 


- H 2* D E L A 


+ 


H 2 


* U A 


D L T RDD = 


-H1*DLTRD 


- H 2* D L T R 


+ 


H 2 


* U T 


T H R 0 T D = 


- UE*T HROT 


+ K*UE*D LTR 







MOMENT EQUATIONS 



FORCE EQUATIONS 



I NTE GRAT I ONS 

P = I NTGRL ( PO , PD ) 

Q = INTGRL(QO.QD) 

R = INTGRL(RO,RD) 

UU= I NTGRL (UO , UD ) 

VV= I NTGR L ( VO , VD ) 

WW= I NTGRL (WO , UD ) 

,******************** EULER ANGLES 
PHI = I NTGRL (PH I O', PH ID) 

THT = INTGRL(THTO,THTD) 

PS I = INTGRL(PSIO,PSID) 

P H I D E G = RAD EG* PH I 
T H T D E G = RADEG*THT 
P S I D E G = RADEG*PS I 

,******************** SERVO EQUATIONS 
DELED = SERVO(DELEDO,DELEDD,RMAX) 
DELRD = SERVO(DELRDO,DELRDD ,RMAX) 
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DELAD = SER VO ( D E L AD 0 , DE LADD , RMAX ) 

DLTRD = SER VO ( D LT RDO , D LT ROD , RMAXT ) 

DELE = I NTGRL (DELEO , DELED ) 

DELR = SERVO( DELRO , DELRD , MAXD F L ) 

DELA = INTGRL(DELAO, DELAD) 

DLTR = SERVO (D LTRO , DLTRD , TAMAX ) 

THROT = I NTGRLCTHROTO , THROTD ) 

*** ALTITUDE 

ALT = INTGRL(ALTO,UERTH) 

D I S T = I NTGRL ( GR NSPO , GR NSPD ) 

* 

CONTRL F I N T I M = 5 

PRINT .5'i*'PHIDEG,THTDEG,PSIDEG,P,Q,R,UU f VV,UU,DELA 
SAVE .05, PHIDEG, T H T D E G , P S I D E G , P , Q , R , UU , V V , UU 
PARAM UE = 0 . , U R = 0 . , U A = • . 0 5 2 , U T = 0.,DELA0 = 0 
LABEL ( G 1 ) TIME RESPONSE TO UA =-.052(NO CONTROL) 

GRAPH (G1 ,DE=TEK618) T I M E ( N I = 5 , U N = S E C ) P H I D E G ( N I = 4 , S C = 4 5 , L I = 1 
UN= 1 THTDEG AND PSIDEG ( D E G R E E S ) 1 , L 0 = - 9 0 ) T H T D E G ( A X = 0 M I T , N I = 4 
SC=45,LI=3,LO=-90) PSIDEG(AX=OMIT,NI=4,SC=45,LI=4,LO=-90) 
GRAPH < G2 , DE = TEK61 8 , OV , P0 = 0 , 5 ) T I M E ( N I = 5 , U N = S E C ) P(NI=4, ... 

U N = 1 Q AND R ( R AD /SE C ) 1 , LO= - 2 , SC= 1 , L I = 1 ) Q(AX = OMIT, ... 
NI=4,LO=-2,SC = 1 , L I = 3 ) R(AX = OMIT,NI=4,LO=-2,SC = 1 , L I = 4 ) 
GRAPH (G2,DE = TEK618,OV, PO = 6, 0) T I M E ( N I = 5 , U N = S E C ) VV(NI=4, ... 

U N = 1 UU AND UU ( FT/SEC ) 1 , LO= - 2 , SC = 5 , L I =3 ) UU(AX = OMIT, .. 
NI=4,LO=-2,SC=5,LI=1 ) UU(AX=OMIT,NI=4,LO=-2,SC=5,LI=4) 

END 

STOP 
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APPENDIX C 



OCNERDIIER AND LINEAR MDDEL SUUIAITCN 

Hie DSI/VS source code used to simulate the linearized model is 
listed here. The steady-state controller is implemented under the 
"SAMPLE" segment. This same controller is used in the nonlinear model 



TITLE SIMPLE MODEL WITH OPTIMAL GAINS FROM CONTROLS - 9 AUG 87 

*** FEED BACK GAINS USING THE ERROR STATE CONTROL LAW U= -F(X-R) *** 

ARRAY RKPT(8),RKQR(8),XPT(8),XQR(8),UPT<2),UQR(2),EPT(8),EQR(8) 

FGEN2 F B K P T = 2,8,16, 1 , 2 , 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , ... 

-1.994 84, -.60299, -.71717, .27437, 1.74559, .00994, .14049, .00016, .. 

-.58535, 111. 16428, -.20533, 28. 90878, .58170, 1.66845, .00962, .13667 

* 

FGEN2 F B KQR = 2,8,16, 1 , 2 , 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , ... 

-1.61288, -1.61027, -.55320, .31985, 1.05412, -.20 75 2, .12 65 6, -.00264,. 
1.5 6829, -1.5 9987, -.304 79, -.5 66 15, .10693, 1.23477, .000 96, .12 962 

★ 

CONST RMAXT = 100., T AM A X = 100., MAXD F L = . 5236, R M A X = .87266 

CONST A1 = . 0236133, A2 = -3. 7381, A3=-6. 750686, A4=. 4830933 
CONST A5_=_6 . 784433, A6=. 107891, A7=. 893727 
CONST LAE = -21 .29,MEE = - 14.51 ,NRE = - 16.68 
CONST H 1 = 17.77, H2= 157.91 

CONST PHDES=.0,ALRDES=.00,THDES=.0,PSDES=.0 

CONST RKPT ( 3 ) = 0 , RKPT (4)=0, RKPT (5 ) = 0 , RKPT (6) = 0 , RKPT ( 7 ) = 0 , RKPT (8)=0 
CONST RKQR(3)=0, RKQR<4)=0, RKQR<5 )=0, RKQR(6)=0, RKQR(7)=0 , RKQR(8)=0 
I N I T 

D E L E D 0 = 0 
D E L R D 0 = 0 
D E LAD 0= 0 
D L T R D 0 = 0 
D E L E 0 = 0 
D E L R 0 = 0 
D E L A 0 = 0 
D L T R 0 = 0 
P0= 0 
Q 0 = 0 
R 0 = 0 
T H R 0 = 0 
1= 1 
J= 1 
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********* define macro for use in servo equation limiting 

MACRO OU T = SERVOCOUTO , OUTD , L I M) 

OU T = I NTGRL (OUTO , OUTD ) 

OU T = LIMIT(-LIM,LIM,OUT) 

E NDMAC 

DYNAMI C 

**** INPUT CALCULATIONS 

RKPT<1) = PHD E S* ST E P < 0 . ) 

R K P T ( 2 ) = ALRDES*STEP<0 . ) 

RKQR(1) = T HD ES*ST EP < 0 . ) 

RKQR ( 2 ) = PSDES*STEP(0. ) 



X P T ( 1 ) = PH I B 
X PT ( 2 ) = U 
XPT (3 ) = P 
X P T ( 4 ) = THROT 
XPT ( 5 ) = DELA 
XPT (6) = DLTR 
XPT ( 7 ) = DELAD 
XPT ( 8 ) = DLTRD 
XQR ( 1 )= THTB 
XQR ( 2 )= PS I B 
X Q R ( 3 ) = Q 
X Q R ( 4 ) = R 
XQ R ( 5 ) = DELE 
XQ R ( 6 ) = DELR 
X Q R < 7 ) = DELED 
X Q R ( 8 ) = DELRD 



SAMPLE 



D Q_3 0 1 = 1,2 
BUPT =0 
BUQR =0 
DO 20 J = 1 , 8 

EPTC J )= XPT ( J ) - RKPT ( J ) 

EQR(J)= XQR(J) - RKQR(J) 
F B P T = FGEN2 ( FBKPT , I , J ) 

F B Q R = FGEN2C FBKQR , I , J ) 

B U P T = BUPT + EPT (J )*FBPT 
B U Q R = BUQR + EQR(J)*FBQR 
20 CONTINUE 

UPT ( I ) = • (BUPT ) 

UQR ( I ) = - < BUQR ) 

30 CONTINUE 



U A = U P T ( 1 ) 
U T = UPT ( 2 ) 
UE = UQR< 1 ) 
UR = UQR ( 2 ) 
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DERIV 



PH I BD= P 
UD = A7MHR0T 
THTBD= Q 
P S I B D = R 

P D = A 1 * U + A2*THR0T + L A E *D E L A 
T H R 0 T D = -THROT + . 5 * D L T R 

QD = A3 * R + MEE*DELE 
R D = A5*Q + NRE*DELR 
ELEVATOR, RUDDER, AILERON, THROTTLE 



D E L A D D = 


• H 1 * D E L AD 


- H 2 * D E L A 


+ 


H 2 


* U A 


D LTRDD = 


- Hi *DLTRD 


- H2*DLTR 


+ 


H 2 


* U T 


D E L E D D = 


- H 1 *DELED 


- H 2 * D E L E 


+ 


H 2 


*UE 


D E L R D D = 


• H 1 *DELRD 


• H 2 * D E L R 


+ 


H 2 


* U R 



P H I B = INTGRL(PHIBO,PHIBD) 

U =INTGRL(UO,UD) 

T H T B = I NTGRL ( THTBO , T HTBD ) 

P S I B = I NTGRL < PS I BO , PS I BD ) 
P= I NTGRL ( PO , PD ) 

T H ROT = I NTGRL (THRO , THROTD ) 
Q= INTGRL (QO , QD ) 

R = INTGRL (RO , RD ) 



DELED = SERVO(DELEDO , DELEDD , RMAX ) 

DELRD = SE RVO ( D EL RD 0 , D EL RDD , RM AX ) 

DELAD = SERVO(DEL ADO , DE LADD , RMAX ) 

DLTRD = SERVO(DLT RDO , D LTRDD , RMAXT ) 

DELE = INTGRL(DELEO, DELED) 

DELR = SERVO(DELRO, DELRD, MAXDFL) 

DELA = I NT GRL ( D E L AO , D EL AD ) 

QULR = SERVOCDLTRO, DLTRD, TAMAX) 

D E S = . 1 
Z E R = 0. 

CONTROL FINTIM=5,DELS=.04,DELT=.005 
PARAM PHIB0=.00,U0=0,THTB0=0,PSIB0=0 
SAVE .04,PHIB,U,THTB,PSIB,ZER,DES 
PRINT .2,PHIB,U,THTB,PSIB 

LABEL ( G 1 ) REGULATOR GAINS - ERROR STATE - PHI <0)= .1 

GRAPH < G 1 ,DE = TEK618) T I M E ( N I = 5 , U N = S E C ) P H I B ( N I = 4 , S C = . 0 5 , L I = 1 , ... 

UN = 1 R/SEC AND ALT RATE F / S E C ' , L 0 = ■ . 0 5 ) U < A X = 0 M I T , N I = 4 , ... 

SC= . 5 , L I =3 , L0= - . 5 ) DES< AX = OM I T , N I =1 0 , SC = . 05 , L I =4 , L0= - . 05 ) ... 

ZER(AX=OMIT,NI=4,SC=.05,LI=4,LO=-.05) 

GRAPH ( G2 , DE = TEK61 8 , OV , P0 = 6 , 4 . 5 ) T I M E ( N I = 5 , U N = S E C ) T H T B ( N I =4 , . . . 

U N = 1 A N D PSIB R/SEC ' , L0= - . 05 , SC= . 05 , L I = 1 ) P S I B ( AX = OM I T , ... 

N I =4 , LO= - .05,SC=.05,LI=3) DES(AX = OMIT,NI=4,LO=-.05,SC=.05,LI=4) 
ZER(AX=OMIT,NI=4,SC=.05,LI=4,L0=-.05) 

END 

STOP 
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APFQOIX D 



NUMERICAL EUX3K3-BQJETICS OOCE 



1. Modelling the ARPD Body 

The Numerical Electromagnetics Code (NEC) enables the user to 
simulate a solid surface using either the wire grid or the surface patch 
technique. Extensive use of symmetry in the program and the circular 
symmetry of the AROD model allcw that the surface patch approach be 
easily implemented. Figure D.l illustrates the procedure schematically. 




Figure D.l Schematic for AROD Body Modelling 
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Hie following computer code is the NEC model of the body which was 
utilized. 

(1) Build the 1st column of patches. The column is 9.92° in angular 
width (arc) and 14 inches high. The arc of the column is in the 
XY plane and the height of the column is in the Z direction at a 
radius of 15 inches from the Z axis. 



SP 


0, 


0, 


15. , 


o. , 


5. , 


1 8 . 435 , 


0. , 


9 . 92 


SP 


0, 


0, 


15. , 


o. , 


1.5, 


5.711, 


0. , 


7.44 


SP 


0, 


0, 


15. , 


o. , 


•1.5, 


-5.711, 


o. , 


7.44 


sp" 


0, 


0, 


15. , 


o. , 


•5. , - 


1 8 . 435 , 


o. , 


9.92 



(2) Rotate and duplicate the column. The 1st column is rotated in the 
XY plane and duplicated to form a column that is now 18.84° wide 
and 14 inches high. 

GM 0, 1, 0., 0., 9.4735, 0., 0., 0. 

(3) Add patches to top and bottom to form wedge. These patches 
together with the column form a 18.84° section of the cylinder. 



SP 


0, 


0, 


13.95, 


1.16, 


7 . , 


90. , 


4.737, 


9.26 


SP 


0. 


0, 


13.95, 


1.16, 


- 7. , 


-90. , 


4 . 737, 


9.26 


SP 


0, 


0, 


11.96, 


,990, 


7 . , 


90. , 


4.737, 


7.94 


SP 


0. 


0, 


11.96, 


. 990 , 


-7. , 


-90 . , 


4 . 737, 


7.94 


SP 


0, 


0, 


9.470 , 


. 785 , 


7. , 


90. , 


4 . 737, 


9.425 


SP 


0, 


0, 


9.470, 


. 785 , 


-7. , 


-90. , 


4 . 737, 


9.425 


SP 


0, 


0, 


4.980, 


.413, 


7. , 


90. , 


4.737, 


9.92 


SP 


0. 


0. 


4 . 980 , 


.413, 


- 7 . , 


- 90 . , 


4 . 737, 


9.92 



(4) Complete the cylinder. The entire 18.84° wedge is rotated 19 
times to form a 360° cylinder which is missing only a top and 
bottom center piece. These two center pieces are added along with 
each antenna that is to be simulated, gs scales the dimensions to 
meters from inches and ge declares an end to geometry. 

GR 0, 19 
GS 2 
GE 

(5) Specify frequency. The frequency at which the model and antenna 
are to simulated is specified. In this example, f = 800 MHz. 

FR 0,0,0,0,800. ,o. 

(6) Store the model. The entire symmetric model is stored in a 
Green's function matrix for recall when the antenna is simulated. 

WG 
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2. Antenna Simulation 



A dipole array and a loop antenna are simulated in the NEC code 
which follows. 

a. 5 Dipole Ring Array 

A five dipole ring array consists of five half-wave dipoles that 
are arranged in a circle. 

(1) Read the Green's function. The AROD body model which is stored in 
a Green's function matrix is recalled to start the simulation. 

This matrix also includes the frequency at which the antenna will 
transmit. 

GF 

(2) Specify the dipoles. Wire .1 inch in radius and 7.38 inches long 
is used to model a half-wave dipole at f = 800 MHz. 11 segments 
in each dipole is used so that NEC automatically calculates 
current losses. 5 dipoles are spaced equally around the circle at 
a radius of 15.5 inches. The dipoles longitudinal axis is in the 
Z direction. 



GU 


i , 


11 , 


15.5, 0., 3.69, 15.5, 0.,-3.69, 


. 1 


GW 


2, 


11 , 


4.79, 14.74, 3.69, 4.79, 14.74,- 


3.69, .1 


GU 


3, 


11 , 


-12.54, 9.11, 3.69, -12.54, 9.11 


,-3.69, . 1 


GU 


A, 


11 . 


-12.54, -9.11, 3.69, -12.54, -9. 


11,-3.69, 


GU 


5, 


11 , 


4.79,-14.74, 3.69, 4.79, -14.74, 


-3.69, .1 



(3) Complete the cylinder. The two center patches which were excluded 
from the symmetric AROD body model are added, the dimensions are 
scaled to meters, and the geometry is ended with no ground plane. 

SP 0 , 0 , 0 . , 0 ... 7 . , 90 , 0 . , 12.57 
SP 0, 0, 0. , 0. , -7. , -90, 0. , 12.57 
GS 2 
GE 0 

(4) Specify the driving source. In this example, a feed point at the 
center of each dipole provides 1.0 volts. 

ex o , i , 6 , oi , i . , o . 



b. loop Antenna 

A single wire loop antenna is added to the AROD body model. The 
loop is flush to the surface of the AROD. 
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(1) Read the Green's function. The AROD body model which is stored in 
a Green's function matrix is recalled to start the simulation. 

This matrix also includes the frequency at which the antenna will 
transmit. 

GF 

(2) Specify 1st wire in the circle. The circle will consist of wire 
arc segments connected together. The 1st segment spans 10° and is 
at a radius of 15.5 inches. The wire is in the XZ plane initially 
and is then moved to the XY plane. 

GA 1 , 3 , 15.5, - 5 . , 5 . , .1 

GM 0, 0, 90., 0., 0., 0., 0., 0., 1 

(3) Create the loop. The first wire is duplicated and rotated 36 
times to form a complete circle. 

GR 1, 36 

(3) Complete the cylinder. The two center patches which were excluded 
from the symmetric AROD body model are added, the dimensions are 
scaled to meters, and the geometry is ended with no ground plane. 

SP 0 , 0 , 0 . , 0 . , 7. , 90 , 0 . , 12.57 

SP 0, 0, 0., 0 . , * 7 . , - 90 , 0., 12.57 

GS 2 
GE 0 

(4) Specify the driving source. In this example, a feed point at (p = 
0° provides 1.0 volts. 

ex o, 1, 6, 01, i., o. 
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